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Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement. 


To  Convert  From 


angstrom 

atmosphere  (normal) 

bar 

bam 

British  thermal  unit  (thermochemical) 
calorie  (thermochemical) 

2 

cal  (thermocheinical)/cm 
curie 

degree  (angle) 
degree  Fahrenheit 


meters  (m) 
kilo  pascal  (kPa) 

kilo  pascal  (kPa) 

2 2 
meter  (m  ) 

joule  (J) 

joule  (J) 

2 2 
mega  joule /m  (MJ/m  ) 

•giga  becquerel  (GBq) 

radian  (rad) 


Multiply  By 


1.000  000  X E -10 
1.013  25  X E 42 
1.000  000  X E 4 2 
1.  000  000  X E -28 
1.  054  350  X E 43 
4. 184  000 
4.  184  000  X E -2 
3.  700  000  X E 4l 
1.  745  329  X E -2 
t = (t*f  4 459.67)/!.  8 


electron  volt 

joule  (J) 

1.602  19  X E -19 

erg 

joule  (J) 

1 . 000  000  X E -7 

erg /second 

watt  (W) 

1 . 000  000  X E -7 

foot 

meter  (m) 

3.  048  000  X E -1 

foot  -pound  -force 
gallon  (l*  S.  liquid) 

joule  (J) 

, 3.3. 

meter  (m  ) 

1 1.355  818 

3.  785  412  X E -3 

inch 

meter  (m) 

2 54  0 000  X E -2 

jerk 

joule  (J) 

1 000  000  X E * 9 

j joule Ailogram  (J /kg)  (radiation  dose 
absorbed) 

Gray  (Gy) 

1 . 000  000 

kilotons 

terajoules 

4.  183 

kip  (1000  Ibf) 

newton  (N) 

. 4.  448  222  X E 43 

kip /inch2  (ksi) 

kilo  pascal  (kPa) 

6 894  757  X E 43 

ktap 

newton -second  /m 
(N-8  /m2) 

1 . 000  000  X E 4 2 

micron 

meter  (m) 

1 000  000  X E -6 

mil 

meter  (m) 

2.  540  000  X E -5 

mile  (international) 

meter  (m) 

1. 609  344  X E 43 

ounce 

kilogram  (kg) 

2. 834  952  X E -2 

pound -force  (lbs  avoirdupois) 

newton  (N) 

4. 448  222 

pound -force  inch 

newton -meter  (N*m) 

1.  129  848  X E -1 

pound  -force /inch 

newton /meter  (N/m) 

1.  751  268  X E 42 

pound  -force  Aoot 

kilo  pascal  (kPa) 

4.  788  026  X E -2 

pound -force /incli  (psi) 

kilo  pascal  (kPa) 

6.  894  757 

pound -mass  (lbm  avoirdupois) 

kilogram  (kg) 

4.  535  924  X E -1 

pound -mass -foot  (moment  of  inertia) 

kilogram -meter 
(kg-m2) 

4.  214  011  X E -2 

| pound -mass /Foot3 

kilogram  /mete  r‘ 

(kg/m3) 

1.601  846  X E 4l 

rad  (radiation  dose  absorbed) 

••Gray  (Gy) 

1.  000  000  X E -2 

roentgen 

| 

coulomb  Ai  logram 
(C/kg) 

2.  579  760  X E -4 

shake 

second  (s) 

1 . 000  000  X E -8 

slug 

kilogram  (kg) 

1.459  390  X E 4l 

torr  (mm  Hg,  0*  C) 

kilo  pascal  (kPa) 

1.333  22  X E -l 

•The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1 Bq  - 1 event/s. 

••The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 

A more  complete  listing  of  conversions  may  be  found  in  "Metric  Practice  ( 
American  Society  for  Testing  and  Materials. 
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SECTION  I 
INTRODUCTION 


This  study  is  concerned  with  the  simulation  of  the  situation  where 
an  aircraft  inlet  is  subjected  to  a blast  wave  of  sufficient  strength 
and  long  enough  duration  to  produce  significant  distortion  of  the  flow 
field  at  the  engine  face.  For  critically  designed  engine- inlet  systems, 
such  blasts,  produced  by  nuclear  explosions,  might  cause  severe  engine 
stall  or  flameout  problems. 

The  objective  of  the  test  program  reported  here  was  to  develop  a 
test  technique  for  measuring  blast  interaction  with  an  aircraft  inlet 
and  to  conduct  a series  of  tests  for  the  development  of  prediction 
methods . 

Wind  tunnels  provide  a well-developed  means  for  the  simulation  of 
flight  flow  conditions  for  an  engine  inlet.  A technique  for  the  simulation 
of  blast  waves  in  a supersonic  wind  tunnel  by  employing  a shock  tube 
mounted  in  the  wind  tunnel  wall  has  been  demonstrated  by  Pierce  (Ref.  1.1). 

A development  and  test  program,  conducted  at  the  Arnold  Engineering 
Development  Center  (Ref.  1.2),  is  reported  here  for  the  simulation  of 
blast  wave  intercepts  with  a scaled  aircraft  inlet  in  subsonic  flight 
using  the  shock-tube  technique  for  firing  the  blast-type  waves.  Pilot 
tests  of  the  technique  were  performed  in  a 1-ft  wind  tunnel  during 
August  1975,  followed  by  "bench"  tests  of  a prototype  shock  tube  during 
August  1976.  A wind  tunnel  test  program  then  was  conducted  in  the  AEDC 
16T  wind  tunnel  with  a 0.1-scale  B-l  inlet  during  September  24-28, 

1976. 

This  report  is  concerned  primarily  with  documentation  of  the  development 
of  the  facility  to  test  the  response  of  engine-inlet  systems  to  blast 
waves  in  high  speed  wind  tunnels,  with  the  presentation  of  test  results 
from  the  facility  for  a 0.1-scale  aircraft  inlet,  and  with  the  discussion 
and  application  of  associated  analytical  techniques  to  predict  inlet 
behavior  in  a blast  environment. 

Section  2 of  the  report  describes  preliminary  analytical  studies 
and  small  model  wind  tunnel  tests  performed  to  assess  the  feasibility  of 
performing  well  controlled  inlet  blast  tests  in  a high  speed  wind  tunnel 
with  blast  waves  produced  by  shock  tubes  fired  into  the  tunnel.  Section 


3 describes  the  subsequent  development  of  a large  shock  tube  facility 
for  use  in  the  AEDC  16T  (16  foot  square)  transonic  wind  tunnel  (Ref. 

1.2).  Section  4 describes  and  presents  sample  results  of  a blast  test 
program  performed  in  this  facility,  using  three  23  in. -diameter  shock 
tubes  and  a 0.1-scale  B-l  inlet  model.  Section  5 presents  a general 
analysis  of  the  test  results,  pointing  out  some  basic  characteristics  of 
the  inlet  blast  interaction  process  and  presenting  an  analysis  of  the 
inlet  distortion  measurements  obtained  during  the  tests. 

Section  6 presents  a comparison  of  experimental  blast  pressure  time 
histories  with  theoretical  calculations  made  with  the  BID  computer  code. 
Section  7 supplements  the  test  results  by  presenting  some  blast  calculations 
made  with  the  BID  code  to  illustrate  some  blast-inlet  interaction  effects 
over  a greater  range  of  blast  intercept  angles  than  were  covered  by  the 
tests. 

Section  8 discusses  the  interaction  of  a blast-induced  inlet  shock 
with  an  engine  fan  and  Section  9 discusses  blast-induced  boundary  layer 
interactions  in  an  inlet. 

The  overall  study  results  are  discussed  in  Section  10  and  conclusions 
are  given  in  Section  11. 

Appendix  A discusses  some  theoretical  studies  of  the  firing  of  a 
shock  tube  into  a stationary  fluid.  Appendix  B presents  the  steady- 
state  pre-blast  measurements  made  during  the  16T  tests. 


SECTION  II 

PRELIMINARY  INVESTIGATION  FOR  DEVELOPMENT  OF  TEST  TECHNIQUE 


As  the  first  stage  of  the  development  of  a wind-tunnel  shock-tube 
blast  test  facility,  a program  of  preliminary  theoretical  calculations 
and  model  tests  was  undertaken  to  establish  the  feasibility  of  the  test 
concept  and  to  provide  design  information  for  construction  of  a full- 
scale  installation. 


2.1  THEORETICAL  STUDIES-REFLECT/S2D 

The  first  development  step  consisted  of  theoretical  calculations 
for  the  problem  of  a circular  shock  tube  fired  into  stationary  air. 

These  calculations,  made  with  the  Kaman  AviDyne  REFLECT-2  and  S2D  computer 
codes  (see  Appendix  A)  indicated  that  useful  shock  wave  patterns  simulating 
nuclear  blast  wave  patterns  in  the  2 to  5 psi  overpressure  range  could 
be  obtained  at  distances  3 to  10  diameters  downstream  of  the  end  of  a 
shock  tube.  In  particular,  these  calculation  results,  summarized  in 
Appendix  A,  indicated  that,  at  a distance  of  about  4.5  diameters  down- 
stream of  the  end  of  a shock  tube,  a blast  duration  of  nearly  constant 
overpressure  could  be  obtained  for  about  2.5  milliseconds  per  foot  of 
tube  diameter. 

2.2  SHOCK  TUBE  FIRINGS  IN  AEDC  IT  WIND  TUNNEL 

Using  the  results  of  the  above  theoretical  calculations  as  a 
guideline,  a small  2 in.-dia  shock  tube  was  constructed  and  was  fired 
into  the  AEDC  IT  (one-foot-square)  transonic  wind  tunnel  for  a wide 
range  of  flow  and  shock  tube  conditions,  during  which  pressure  measurements 
were  made  for  a large  number  of  locations  inside  the  wind  tunnel. 

The  geometric  setup  used  in  the  AEDC  IT  tunnel  is  indicated  in 
Figure  2.1.  Three  shock  tube  configurations  were  tested,  about  14 
inches  long,  with  ratios  of  the  driven-tube  lengths  to  driver-tube 
lengths  of  about  0.5,  2 and  8.  The  end  of  each  shock  tube  protuded  into 
the  tunnel  by  about  one  inch  (perpendicular  to  the  floor)  where  it  was 
terminated  by  a 8 in-dia.  circular  flange  (Fig.  2.1).  Firings  were  made 
with  the  angle  of  the  shock  tube  axis  to  the  tunnel  axis,  $,  (see 
Figure  2.1),  being  90°,  60°  and  45°.  Tunnel  ambient  pressures  were  in 
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Figure  2.1.  Test  setup  for  blast  tests  in  the  AEDC  IT  transonic  wind  tunnel  (tunnel  side  walls  removed). 
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the  range  11  to  16  psia  and  shock  tube  driver  pressures  ranged  from 
about  60  to  380  psia. 

Transient  blast  pressure  measurements  were  made  inside  the  wind 
tunnel  by  means  of  pressure  transducers  imbedded  in  the  surface  of  a 
thin  plate  placed  in  the  plane  of  symmetry,  as  shown  in  Figure  2.1. 
Measurements  were  made  of  local  blast  pressure  time  histories  for  radial 
distances  from  the  shock  tube  exit  of  2.5,  3.5  and  4.5  tube  diameters 
(r/d)  and  for  polar  angle  locations  (0  in  Fig.  2.1)  of  -15°,  0°,  15°, 

30°,  45°,  60°  and  75°.  Also,  for  several  locations,  measurements  of 
transient  total  pressure  and  flow  velocities  were  made  with  two  total- 
pressure  and  two  claw  probes  mounted  about  a half  inch  from  the  face  of 
the  thin  plate,  as  shown  in  the  figure.  Kulite  XCQL-093-025  pressure 
transducers  were  used  in  the  plate  and  probes. 

Sample  shadowgrams  showing  the  blast  field  produced  by  firing  the 
2-in  shock  tube  into  the  wind  tunnel  are  presented  in  Figures  2.2  and 
2.3.  Figure  2.2  shows  three  shadowgrams  of  one  firing  into  still  air  and 
Figure  2.3  shows  two  shadowgrams  for  a firing  into  a Mach  0.8  tunnel  flow. 
The  shadowgrams  are  taken  looking  into  the  side  of  the  tunnel,  with  the 
tunnel  flow  proceeding  from  left  to  right,  and  with  the  shock  tube  being 
fired  into  the  tunnel  vertically  up  (<J>=90°)  through  the  bottom  wall 
of  the  tunnel. 

In  Figure  2.2,  the  lowest  shadowgram  was  taken  first  (before  blast 
exit  from  the  tubes)  and  the  upper  one  last.  The  blast  front  can  be 
clearly  seen  in  these  shadowgrams  and  becomes  essentially  a circular 
arc  after  the  front  of  the  blast  wave  has  moved  several  tube  diameters 
away  from  the  tube  exit.  In  the  upper  shadowgram  of  Figure  2.2,  between 
the  top  of  the  blast  front  and  the  shock-tube  exit  is  seen  a nearly 
horizontal  dark  line  which  roughly  delineates  the  contact  surface  between 
the  "hot"  tunnel  air  (above)  and  the  "cold"  expanded  air  (below)  from 
the  driver  section  of  the  shock  tube.  This  cold  gas  can  also  be  more 
clearly  seen  as  the  dark  jet  core  in  Figure  2.3.  From  the  viewpoint  of 
simulating  a nuclear  blast  or  other  long-duration  blast  wave,  it  should 
be  noted  that  only  that  portion  of  the  blast  wave  lying  between  the  blast 
front  and  the  hot-cold  gas  contact  surface  is  of  interest. 


15 


The  effect  of  the  tunnel  flow  on  the  shape  of  the  blast  shock  front 
was  found  to  be  relatively  small  in  the  sense  that  the  blast  front  shape 
remained  essentially  circular  for  all  Mach  numbers  tested  (up  to  0.9), 
but,  of  course,  the  pattern  is  convected  downstream  more  with  increasing 
tunnel  speed.  The  central  cold-gas  core  also  tends  to  be  convected 
downstream  (Fig.  2.3),  but  less  so  than  the  blast  front. 

Pressure  time  histories  obtained  from  all  transducers  were  digitized 
by  ARO  and  were  subsequently  provided  in  the  form  of  plots,  tabulations 
and  magnetic  tapes  (Ref.  2.1). 

Sample  pressure  measurements  in  the  blast  field  are  shown  in 
Figure  2.4  for  the  shock  tube  normal  to  the  tunnel  axis,  $=90°,  measured 
at  a radial  distance  of  4.5  tube  diameters  from  the  tube  exit,  at  a 
tunnel  Mach  number  of  0.8  and  a ratio  of  tube  driver  pressure  to  driven 
pressure  of  about  10.  It  should  be  noted  that  for  a polar  angle  (0 
in  Figure  2.1)  of  15°  a blast  wave  is  produced  with  a nearly  constant 
pressure  level  for  about  0.2  milliseconds.  In  the  polar  angle  range  - 
15°  to  30°  useful  nearly-constant  blast  durations  of  at  least  0.1  millisecond 
are  seen  to  have  been  produced;  but  for  larger  angles  (0  > 45°)  the 
pressure  has  an  initial  spike-behavior  which  makes  it  of  little  value 
for  simulating  blast  shocks  of  long  duration. 

Similar  results  to  those  shown  in  Figure  2.4  were  obtained  for  most 
conditions  tested,  covering  shock-tube  angles  (4>)  of  45,  60  and  90 
degrees  to  the  tunnel  axis,  tunnel  Mach  numbers  of  0,  0.6,  0.75  and  0.9 
and  driver/driven-tube  pressure  ratios  between  about  4 and  30.  Useful 
nearly-constant  pressure  blast  durations  up  to  a maximum  of  about  0.3 
millisecond  were  obtained  for  some  conditions. 

In  order  to  provide  a more  quantitative  index  of  the  obtained 
duration  of  nearly  constant  pressure,  estimates  were  made  for  the  90° 
tube  angle  (<J>  = 90°)  from  curves  like  Figure  2.4  of  the  time,  t^o’ 
required  for  the  blast  overpressure  to  decay  (or  increase)  30  percent 
from  the  initial  shock  value.  For  blast  overpressures  between  2 and  5 
psi  at  a radial  distance  of  4.5  diameters  (r/d  = 4.5),  the  30-percent- 
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decay  time  was  found  to  generally  increase  about  linearly  with  increasing 
overpressure.  In  the  4-5  psi  overpressure  range  the  decay  time  was 
about  0.16  msec  or  1 msec/ft.  dia.  for  a polar  angle  (0)  of  0°  and  about 
0.09  msec  or  0.5  msec/ft  for  30°,  and  did  not  appear  to  vary  much  with 
Mach  number.  For  a polar  angle  of  15°  the  decay  time  decreased  considerably 
with  increasing  Mach  number;  at  the  5 psi  level,  it  appeared  to  vary 
from  about  0.3  msec  or  2 msec/ft  at  Mach  0.6  to  about  0.18  msec  or  1.1 
msec/ft  at  Mach  0.9. 

The  angle  of  the  tubes  to  the  tunnel  axis  (<5)  had  a marked  effect 
on  the  blast  waves  that  were  produced.  As  the  angle  was  reduced,  causing 
the  blast  to  be  fired  more  in  the  downstream  direction,  the  resulting 
shock  overpressure  was  lower  and  the  overpressure  at  a fixed  point 
decayed  more  rapidly.  Satisfactory  results  were  obtained  for  the  tube 
at  90  degrees  to  the  tunnel  axis,  so  this  angle  was  selected  for  later 
development  of  a large  shock  tube  facility. 

Blast  pressure  waveforms  were  found  to  be  relatively  insensitive  to 
the  shock-tube  driven-tube/driver  length  ratio  in  the  useful  early-time 
(hot-gas)  part  of  the  flow  period.  Hardly  any  difference  was  observed 
between  the  waveforms  for  the  length  ratios  of  0.5  and  2.0,  but  some 
waveforms  for  the  length  ratio  of  8 had  somewhat  shorter  nearly-constant 
durations  than  those  for  the  other  length  ratios. 

To  summarize  the  test  results,  it  appeared  that  useful  long  duration 
blast  waves  could  be  obtained  by  firing  shock  tubes  into  a wind  tunnel 
operating  at  speeds  up  to  at  least  Mach  0.9.  Nearly  constant  pressure 
blast  durations  lasting  up  to  a maximum  of  about  2 milliseconds  per  foot 
of  tube  diameter  were  obtained  for  some  firing  conditions  and  the  longest 
duration  of  nearly  constant  pressure  blast  waves  was  obtained  at  polar 
angles  (0)  from  the  shock  tube  of  about  15°. 
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SECTION  III 

DEVELOPMENT  OF  16T  BLAST  TEST  FACILITY 

Having  established  the  feasibility  of  shock  tube-wind  tunnel  blast 
tests  by  small  model  tests  in  the  IT  tunnel,  a program  was  undertaken  to 
develop  and  apply  a similar  large  shock  tube  facility  for  engine  inlet 
testing  in  the  AEDC  16T  (sixteen  foot  square)  transonic  wind  tunnel 
(Ref.  2.1).  Development  of  this  facility  consisted  of  selection  and 
modification  of  a inlet  model,  shock  tube  design,  construction  and 
testing,  selection  of  dynamic  instrumentation  for  measuring  the  blast 
field,  and  tunnel  modifications. 

The  test  installation  was  designed  to  take  as  much  advantage  as 
possible  of  existing  hardware  and  fixtures. 

3.1  GENERAL  DESIGN  CRITERIA 

The  inlet  model  selected  for  blast  testing  in  the  AEDC  16T  tunnel  is 
designated  as  the  0.1-scale  B-l  Inlet  Development  Model  II.  This  model 
was  selected  as  having  a representative  transonic  inlet  of  modern  design 
which  was  already  well  instrumented  with  dynamic  pressure  engine  face 
transducers  of  high  enough  frequencies  to  respond  to  blast  induced 
transients.  The  model  required  only  small  modifications  to  strengthen 
it  for  blast  testing  and  it  was  compatible  with  existing  mounting  fixtures 
for  the  16T  tunnel. 

The  important  features  of  the  model  with  respect  to  shock  tube  design 
are  the  size  of  the  inlet  opening,  which  is  about  5.5"  long  by  4.4"  wide 
by  3.3"  high  for  each  inlet,  and  the  length  of  the  inlet,  about  25". 

Other  model  characteristics  are  discussed  in  Section  4.2. 

The  shock  tubes  had  to  be  designed  to  produce  a sufficiently  small 
variation  of  incident  shock  pressure  across  the  above  inlet  opening. 

Also  the  duration  of  the  blast  wave  produced  had  to  be  at  least  long 
enough  so  that  the  blast  wave  does  not  decay  significantly  before  its 
front  has  had  time  enough  to  pass  down  the  inlet  into  the  engine  and  to 
be  reflected  back  past  the  engine  face. 
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It  appeared  that  the  most  cost-effective  facility  design  satisfying 
the  above  requirements  was  to  position  the  test  model  near  the  center  of 
the  tunnel,  permitting  use  of  existing  fixtures,  and  to  mount  shock 
tubes  in  both  tunnel  walls,  two  in  one  wall  and  one  in  the  other,  as 
indicated  in  Figures  3.1  and  3.2.  This  design  concept,  using  three 
shock  tubes,  permitted  the  firing  of  three  blast  waves  at  the  model 
before  the  tunnel  had  to  be  shut  down  to  change  shock  tube  diaphragms, 
thus  saving  considerable  time  and  expense.  Dimensions  and  orientations 
of  the  three  shock  tubes  were  chosen  to  conform  as  close  as  practical  to 
the  optimum  conditions  obtained  from  the  IT  model  tests,  namely,  ratio 
of  tube-to-inlet  distance  to  tube  diameter  of  about  4.5,  polar  angles 
from  shock  tube  to  inlet  (0)  between  -15°  and  +30°,  with  about  15°  being 
optimum,  and  shock  tube  axis  perpendicular  to  tunnel  wall  (<J>=90°) . A 
common  inside  diameter  of  22.6  inches  appeared  to  be  the  best  choice, 
which  gave  distance/diameter  ratios  of  4.3,  3.8  and  4.3  for  shock  tubes 
1,  2,  and  3,  respectively.  The  resulting  geometrical  arrangement  of 
tunnel,  shock-tubes  and  model  is  shown  schematically  in  Figure  3.1 
and  photographs  of  the  setup  are  presented  in  Figure  3.2. 

With  respect  to  the  question  as  to  how  uniform  the  blast  pressure 
would  be  across  the  inlet  opening,  the  results  of  the  IT  tests  indicated 
that  for  the  above-selected  model  and  shock  tube  geometry,  the  intensity 
of  the  blast  wave  would  generally  vary  across  one  inlet  not  more  than  + 

2 percent  from  the  mean  value  and  at  worst  not  exceeding  + 4 percent. 

3.2  SHOCK  TUBE  DESIGN  AND  TESTING 

Using  the  guidelines  discussed  in  Section  3.1,  ARO  designed  and 

* 

constructed  three  similar  22.6  in-dia  shock-tube  blast  generators  for 

•k 

A necked  down  shock  tube  driven  section  having  a 13.1- inch  exit  diameter 
was  also  designed  as  an  insert  for  any  of  the  three  shock  tubes.  This 
smaller  diameter  insert  was  constructed  and  several  static  firings  were 
made  with  it  to  explore  the  possibility  of  obtaining  a better  blast 
waveform  at  low  blast  pressures  than  could  be  obtained  with  the  large 
diameter  tubes.  This  insert  was  also  intended  to  provide  more  flexibility 
in  the  range  of  blast  overpressures  which  could  oe  obtained  with  a 
particular  pair  of  diaphrams.  However,  the  few  firings  made  with  this 
insert  produced  no  obviously  improved  waveforms  and  the  time  required 
to  install  and  remove  the  insert  from  a shock  tube  when  installed  in 
the  16T  wind  tunnel  proved  to  be  too  long  to  be  tolerated  for  the 
limited  test  time  available.  Consequently,  no  wind  tunnel  firings  were 
made  with  the  small  diameter  driven  section  insert. 
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firing  shock  waves  into  the  AEDC  16T  tunnel.  Some  details  of  one  of 
these  tubes  can  be  seen  in  Figure  3.3.  The  shock  tube  is  about  17  feet 
long.  The  driven-tube  length  is  about  twice  the  driver  tube  length. 

The  tube  has  a double  diaphragm  construction,  and  uses  pre-scored  pre- 
stressed 24  in-dia  diaphragms.  Diaphragms  are  constructed  either  of 
aluminum,  nickel  steel  or  316  stainless  steel,  depending  on  the  firing 
pressure  level. 

All  three  shock  tubes  had  essentially  the  same  design  dimensions 
except  that  the  ratio  of  the  length  of  the  driven  section  to  that  of  the 
driver  section  had  to  be  varied  somewhat  to  permit  close  mounting  of  two 
of  the  tubes  side  by  side  in  the  16T  tunnel  without  interference. 

Judging  from  the  results  of  the  IT  tests  (Sec.  2.2),  these  differences 
in  lengths  would  not  be  expected  to  produce  any  significant  differences 
between  the  blast  waves  generated  by  the  three  tubes. 

Considerable  effort  was  devoted  to  the  design  and  operational  use 
of  the  diaphragms.  A series  of  17  static  firings  was  run  to  establish 
the  range  of  loading  pressures  and  other  requirements  for  achievement  of 
satisfactory  diaphragm  bursts  in  respect  to  the  blast  wave  produced  and 
avoidance  of  ejecta.  A typical  diaphragm  break  pattern  is  shown  in 
Figure  3.4.  The  diaphragms  generally  broke  cleanly  on  the  pre-scored 
lines,  forming  four  petals,  and  never  produced  any  observable  flying 
metallic  debris. 

Shock-tube  operations  were  automated  as  much  as  possible  in  order 
to  obtain  a maximum  number  of  firings  in  the  few  test  days  available  for 
wind  tunnel  testing.  Semi-automated  diaphragm  removal  and  replacement 
required  less  than  15  minutes  per  diaphragm  pair  and  nearly  fully  automated 
pressurization  and  firing  of  a shock  tube  generally  took  less  than  5 
minutes. 

The  above-mentioned  17  static  firings,  made  in  the  AEDC  16S  wind 
tunnel  facility  at  an  ambient  pressure  of  0.5  atmosphere  at  driver 
pressures  from  62  to  233  psia,  were  also  utilized  to  obtain  measurements 
of  blast  overpressures  and  impact  pressures  produced  by  the  shock  tube 
flow,  both  for  protective  design  of  wind  tunnel  equipment  to  be  in  the 
path  of  the  jet  from  the  shock  tube  and  for  providing  a guide  for  diaphram 
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Figure  3.4.  Typical  diaphragm  rupture  pattern 


selection  to  achieve  desired  levels  of  blast  overpressure  at  the  test 
model  during  wind  tunnel  firings. 


3.3  BLAST  INPUT  INSTRUMENTATION 

In  order  to  provide  measurements  of  blast  wave  intensity  and  direction 
as  the  blast  wave  strikes  the  inlet,  three  similar  claw  probes  were 
designed  to  be  mounted  on  or  near  the  test  model.  Figure  3.5  shows  a 
photograph  of  one  of  these  probes  mounted  above  the  tested  inlet.  Each 
probe  consists  basically  of  two  45°-swept  claw  arms  (90°  apart)  with 
Kulite  series  XCQL-093-025  type  dynamic  pressure  transducers  at  the 
tips,  to  measure  total  pressure  and  flow  direction,  and  of  a flat  plate 
surface  toward  the  rear  of  the  probe  with  one  static  pressure  orifice 
and  one  imbedded  Kulite  LQ  series  wafer  pressure  transducer  to  measure 
the  pre-blast  and  transient  "static"  blast  pressures,  respectively.  The 
three  claw  probes  were  calibrated  in  the  AEDC  16T  tunnel  for  sideslip 
angles  from  about  -50  to  + 50  degrees  at  Mach  numbers  of  0.6,  0.7,  0.85 
and  0.95.  The  calibration  data  indicated  that  the  probes  could  provide 
reasonable  estimates  of  flow  characteristics  for  sideslip  angles  between 
about  -30  and  +30  degrees. 

3.4  TUNNEL  WALL  MODIFICATIONS 


Preliminary  KA  estimates  of  the  forces  which  would  be  imposed  on 
the  16T  wind  tunnel  wall  opposite  a firing  shock  tube  indicated  that  the 
local  wall  pressures  could  greatly  exceed  the  tunnel  wall  design  criterion 
of  3 psi  differential  pressure.  To  obtain  design  information  on  such 
wall  loads,  during  the  shock  tube  calibration  test  (Sec.  3.2),  ARO  fired 
one  of  the  22.6  in-dia  shock  tube  at  a large  plate  instrumented  with  an 
array  of  pressure  transducers.  Using  these  test  results  as  a guideline, 
the  tunnel  walls  were  strengthened  in  critical  areas.  In  the  subsequent 
test  program  in  the  wind  tunnel,  no  damage  was  experienced  to  the  tunnel 
walls  during  any  shock  tube  firing. 
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SECTION  IV 

16T  WIND  TUNNEL  BLAST  TEST  PROGRAM 

Following  development  of  the  16T  wind  tunnel  blast  test  facility,  a 
series  of  45  wind  tunnel  shock  tube  firings  were  performed  with  the  0.1- 
scale  B-l  inlet  model.  This  section  presents  a discussion  of  the 
details  of  the  test  setup,  instrumentation,  and  firing  conditions  and 
presents  and  discusses  some  aspects  of  sample  data  obtained  from  the 
tests . 

4.1  GENERAL  TEST  SETUP 

As  was  mentioned  in  Section  3,  the  test  setup  consisted  of  a 0.1- 
scale  B-l  model  installed  in  an  inverted  position  near  the  center  of  the 
AEDC  16T  transonic  tunnel,  with  three  shock  tubes  projecting  through  the 
tunnel  side  walls  (Fig.  3.1).  Figure  3.2  presents  photographs  of  the 
model  and  shock  tube  setup  inside  the  tunnel.  Various  details  of  model 
construction  are  shown  in  Figures  3.5,  4.1  and  4.2.  During  all  tests 
the  model  was  maintained  at  an  angle  of  attack  of  +3  degrees. 

Blast  waves  can  be  fired  at  the  model  inlet  from  any  of  the  three 
22.6  in-dia  shock  tubes  seen  projecting  into  the  tunnel  about  6 inches 
from  the  side  walls  in  Figure  3.2.  The  shock-tube  axes  are  perpendicular 
to  the  side  walls  (<J>=90°)  and  the  axis  of  each  tube  is  located  at  a fixed 
location  between  1 and  4 feet  upstream  of  the  inlet  opening,  so  that, 
for  all  three  shock  tubes,  the  inlet  is  located  within  the  polar  angle 
range  (0)  12  to  29  degrees  from  the  shock-tube  exit  (Figure  3.1),  the 
angles  being  28.3°,  11.2°  and  13.6°  for  tubes  1,  2 and  3,  respectively. 
These  locations  were  chosen  to  give  blast  waves  of  as  long  duration  as 
practical,  whose  fronts  would  generally  strike  the  inlet  roughly  side-on 
or  slightly  from  the  rearward. 

4.2  MODEL  AND  INSTRUMENTATION  DETAILS 

The  most  important  features  of  the  inlet  model  used  for  the  blast 
test  are  illustrated  in  Figures  4.1  through  4.4  (see  also  Ref.  4.1). 

The  model  is  a 0.1-scale  representation  of  the  forward  part  of  the  B-l 
aircraft  fuselage  with  a left-hand  stub  wing  set  at  the  67.5°-  sweep 
position  and  a left-hand  dual  inlet.  The  normally  variable  geometry 


Figure  4.2.  Closeup  view  of  engine  inlet  and  claw  probes  in  the  AEDC  16T  wind  tunnel. 
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Figure  4.3.  Engine  face  transducer  locations.  (Rockwell  drawing) 
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Figure  4.4. 
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inlet  ramp  system  for  the  model  was  fixed  for  these  tests  in  the  normal 
subsonic  cruise  position.  The  engine  face  is  located  at  Nacelle 
Station  269.5.  The  aircraft  engine  is  simulated  by  an  engine  hub 
structure  and  by  remotely  actuated  hydraulically  operated  flow  control 
vanes  which  produce  prescribed  choked  flow  rates  in  the  simulated 
engine  at  a position  about  4.5  inches  downstream  of  the  engine  face. 

Forty  combination  steady-state  and  dynamic  total-pressure  probes 
are  located  in  each  inlet  at  the  engine  face  section,  located  as  indicated 
in  Figures  4.3  and  4.4  and  in  Table  4.1.  The  dynamic  probes  used  Kulite 
type  CQL-080-25  transducers,  which  were  found,  from  pre-test  shock  tube 
calibrations,  to  generally  have  response  times  to  transient  pulses  of 
about  25  microseconds  or  less.  These  Kulite  transducers  had  a 25-psi 
range,  0.080- in  diameter  and  about  125  kHz  natural  frequency. 

Other  model  instrumentation  consisted  of  54  static-pressure  orifices 
and  10  dynamic  Kulite  transducers  whose  locations  are  described  in 
Reference  4.1  (some  indicated  in  Fig.  4.4).  The  10  Kulite  transducers, 
each  accompanied  by  a static-pressure  orifice,  were  mounted  inside  and 
flush  with  the  inlet  wall  in  the  locations  indicated  in  Figure  4.4  (also 
Fig.  5.4)  and  Table  4.1.  These  transducers  were  used  to  monitor  inlet 
ramp  and  cowl  pressures.  All  but  one  of  these  transducers  were  mounted 
in  the  outboard  inlet. 

Also  mounted  on  the  model  were  three  claw  probes,  described  in 
Section  3.3,  to  measure  the  characteristics  of  the  blast  wave  at  the  model. 
Figures  3.5  and  4.2  indicate  the  locations  of  these  three  probes.  Probe 
2 was  located  just  above  the  forward  tip  of  the  inlet;  Probes  1 and  3 
were  located  just  above  the  wing  surface,  about  9.3  inches  outboard  and 
inboard  of  the  inlet  centerline,  respectively,  with  their  tips  located 
at  the  same  longitudinal  station  as  the  inlet  cowl  tips. 

4.3  16T  WIND  TUNNEL  TESTS 

The  16T  wind-tunnel  blast  test  program  was  designed  to  obtain  data 
on  the  blast  response  of  the  test  model  for  a range  of  overpressures 
from  2 to  5 psi  (scaled  to  1 atmosphere  ambient  pressure)  for  tunnel 
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TABLE  4.1 

LOCATIONS  OF  DYNAMIC  PRESSURE  TRANSDUCERS'1" 


Engine  face  pressures  (see  Figure  4.3) 


TYPE 

RADIUS 

(Inches) 

ORIFICE 

NUMBER 

0=22.5 

67.5° 

112.5° 

157.5° 

202.5° 

247.5° 

292.5° 

337.5° 

Total 

1.025 

801 

806 

B9 

mm 

821 

826 

831 

836 

1.393 

802 

807 

mm 

822 

827 

832 

837 

1.683 

803 

808 

mtlm 

823 

828 

833 

838 

1.930 

804 

809 

814 

824 

829 

834 

839 

f 

2.149 

805 

810 

815 

825 

830 

835 

840 

Add  1000  to  tube  number  for  outboard  inlet;  add  2000  to  tube  number  for 
inboard  inlet. 


Cowl  unsteady-state  pressures  (see  Figure  4.4  or  5.4) 


Transducer  Number 

Type 

Location 

1 Outboard 

Inboard 

k 

N.S. 

Remarks 

■M 

2902 

Static 

8.90 

Buzz  indicator 

mm 

- 

ti 

11.60 

Cowl  static 

1904 

- 

it 

16.50 

Cowl  static 

1905 

ti 

25.45 

Cowl  static 

Ramp  unsteady-state  pressures  outboard  inlet  (see  Figure  4.4  or  5.4) 


Mach  numbers  from  0.55  to  0.90  and  for  inlet  flow  rates  representative 
of  cruise  and  maximum  power  conditions.  Table  4.2  outlines  the  test 
matrix  which  was  followed  to  cover  these  conditions. 


Tunnel  testing  was  performed  over  a three  night  period,  on  24 
September  and  27-28  September  1976.  Testing  consisted  of  calibrations, 
steady-state  runs  and  a series  of  blast  firings  as  described  below.  All 
tests  were  performed  at  below  one  atmosphere  ambient  pressures  (usually 
0.5  atmosphere)  to  minimize  blast  loads  on  the  model  and  equipment, 
subject  to  the  limit  of  a minimum  Reynolds  number  per  foot  of  2.5  X 10^. 

4.3.1  Unyawed  Steady-State  Tests 

Prior  to  the  blast  tests,  a series  of  unyawed  steady-state  runs  were 
made  covering  all  Mach  number  and  inlet  mass  flow  rate  ranges  of  interest. 
These  runs  were  made  partly  as  a calibration  of  the  inlet  model  vanes  to 
control  mass  flow  rate  and  partly  to  verify  that  the  inlet  steady-state 
distortion  properties  were  the  same  as  were  found  in  previous  tests  of 
the  model . 

Steady  state  runs  were  also  made  immediately  before  each  shock  tube 
firing  to  provide  reference  (pre-blast)  pressures  at  all  dynamic  pressure 
transducer  locations.  (The  dynamic  pressure  transducers  measure  only 
changes  of  pressure;  absolute  pressures  are  obtained  by  adding  pre-blast 
(steady-state)  pressures  to  the  dynamic  pressures) . Results  of  these 
pre-blast  steady-state  tests  are  presented  in  Appendix  B. 

4.3.2  Yawed  Steady-State  Tests 

A few  steady-state  runs  were  made  with  the  model  at  various  yawed 
positions  in  the  range  -10  to  +10  degrees,  according  to  the  schedule 
given  in  Table  4.3.  These  runs  were  made  partly  to  calibrate  the  inlet 
flow  controls  under  yawed  conditions  and  partly  to  aid  in  assessing  the 
quasi-steady  late-time  response  of  the  inlet  to  the  side-slip  effects  of 
blast  waves  of  very  long  duration. 

4.3.3  Shock  Tube  Firings 

Forty-five  shock  tube  firings  were  performed  in  the  AEDC  16T  tunnel 
at  tunnel  Mach  numbers  of  0,  0.55,  0.70,  0.85  and  0.90.  Tests  were 
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NOMINAL  TEST  CONDITIONS 


Mach  No. 

Flight  conditions 

Mass  f lpw 
(lb/sec  ) 

* 

Overpressure 

(psi) 

Yaw 

(deg) 

0 

Parked 

0 

4 

0 

0.55 

Cruise 

235 

4 

0 

0.55 

Max . power 

350 

4 

0 

0.70 

Cruise 

300 

2,3,4 

0 

0.70 

Max . power 

350 

2, 3, 4, 5 

0,5 

0.85 

Cruise 

300 

2, 3, 4, 5 

0 

0.85 

Max.  power 

350 

4,5 

0 

0.90 

Max . power 

350 

4 

0 

+Full  scale  values  (divide  by  100  to  obtain  model  values) 
* 

Scaled  to  ambient  pressure  of  one  atmosphere 


o o 


TABLE  4.3 


STEADY-STATE  YAW  CONDITIONS 


Mach  no. 

Flow  rate 
(lb/sec) 

Yaw  angle  range 
(deg) 

0.55 

352 

-10  to  +10 

0.70 

303 

-10  to  +10 

0.70 

352 

-10  to  +10 

0.85 

350 

-10  to  +7.5 

generally  performed  with  the  tunnel  static  pressure  set  at  about  1/2 
atmosphere,  except  for  the  Mach  0.55  runs,  for  which  about  2/3  atmosphere 
was  used.  The  inlet  mass  flow  rate  was  generally  set  to  simulate  both 
cruise  and  maximum  power  flight  conditions  for  each  Mach  number  condition 
according  to  the  schedule  in  Table  4.2.  This  schedule  required  testing  at 
mass  flow  rates  of  about  235,  300,  and  350  lb/sec  (full  scale)  . 

Tunnel  steady-state  pre-blast  conditions  are  tabulated  in  Table  4.4 
and  Appendix  B.  Tunnel  total  temperature  was  569°  (+  1°)  for  all  firings 
with  the  tunnel  turned  on  and  was  534°  R for  the  one  static  firing. 

Firings  were  made  over  a range  of  overpressures  from  about  2 to  6 psi, 
scaled  to  a tunnel  ambient  pressure  of  one  atmosphere. 

Also  presented  in  Table  4.4  are  estimated  values  of  the  blast  over- 
pressure (Ap)  and  the  blast  intercept  angle  (<}>)  as  the  blast  wave  strikes 
the  inlet.  The  source  and  accuracy  of  these  values  is  discussed  sub- 
sequently in  Section  4.6. 

4.4  DATA  REDUCTION 

The  test  data  obtained  from  each  shock  tube  firing  were  digitized 
by  ARO  at  time  intervals  of  10.45  microseconds  for  a total  of  2048  samples 
per  firing.  The  resulting  digital  data  were  provided  to  KA  in  the  form 
of  magnetic  tapes  and  time  history  plots  of  most  of  the  basic  data  and 
derived  quantities  (Ref.  4.2).  Tabular  time  history  printouts  of  a few 
variables  were  also  provided  by  ARO. 

4.5  SAMPLE  TEST  RESULTS 

Pressure  time  histories  for  typical  pressure  transducers  are  pre- 
sented in  Figure  4.5  for  a typical  shock  tube  firing  (Run  8/Part  573)  to 
illustrate  the  major  features  of  the  test  data.  In  general,  in  this 
figure,  the  ordinate  label  designates  the  variable  measured  and  the 

vertical  scale  is  always  either  pressure/p  or  pressure/p  (the  latter 

o 

* 

All  mass  flow  rates  in  this  report  are  scaled  to  full  scale  conditions. 
Model  values  are  1/100  of  these  values. 

** 

It  should  be  noted  that  all  blast  pressure  values  designated  as  Ap  in 
this  report  are  all  scaled  to  a tunnel  ambient  pressure  of  one  atmosphere. 
Pressures  not  specifically  designated  as  Ap,  such  as  shock  tube  driver 
pressures  in  Table  4.4  and  all  pressure  values  in  Appendix  B are 
actual  test  pressures. 
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TABLE  4.4 

16T  WIND  TUNNEL  TEST  CONDITIONS 


I 


FIGURE  4.5. 

SAMPLE  TEST  RESULTS 
FOR  RUN  8 (PART  573),  MACH  0.70, 

FLOW  RATE  « 303  LB/SEC,  TUBE  2,  Ap  = 5.0  PSI 


1990  1980  1970 


CR3  PS3 


being  designated  by  an  asterisk  after  the  label),  where  is  the 
steady-state  (pre-blast)  tunnel  total  pressure  and  pQ  is  tRe  steady- 
state  tunnel  ambient  static  pressure. 

The  lowest  part  of  Figure  4.5a  shows  the  blast  pressure  inside  the 
firing  shock  tube  14  inches  from  the  muzzle  end  (STE) . The  pressure 
variation  consists  of  a initial  shock,  followed  by  a slight  rise  and  a 
nearly  constant  level  for  about  8 milliseconds,  after  which  a slow  decay 
follows. 

The  upper  parts  of  Figure  4.5a  give  the  (static)  local  blast  pressures 
(PS)  measured  on  the  three  claw  probes,  which  provide  an  indication  of 
the  strength  of  the  blast  wave  as  it  strikes  the  inlets.  Looking  first 
at  PS2,  which  is  near  to  the  forward  tip  of  the  inlet,  it  may  be  noted 
that  the  pressure  rises  immediately  on  blast  arrival  to  a large  shock 
value,  appears  to  increase  slightly  more  over  a period  of  about  one 
millisecond,  decays  slowly  until  a time  typically  of  about  3.3  milli- 
seconds after  initial  blast  arrival  and  then  decays  more  rapidly. 

Similar  variations  are  seen  for  the  other  two  probes  (PS1  and  PS3) , 
except  that  the  signal  for  PS1,  on  the  blastward  side  of  the  inlet,  is 
complicated  by  the  appearance  of  the  reflection  of  the  shock  wave  from 
the  blastward  blunt  side  of  the  inlet  at  about  one  millisecond  after 
blast  arrival. 

It  should  be  noted  that  for  nuclear  blast  wave  simulation  purposes 
only  the  part  of  the  signal  from  shock  arrival  to  the  rapid  pressure 
change,  starting  at  about  3.3  msec  after  blast  arrival  for  PS2  in 
Figure  4.5a,  is  directly  of  interest.  This  rapid  change  can  appear 
either  as  a rapid  pressure  drop  usually  at  about  3.3  msec.,  as  for  PS2, 
or  it  may  take  the  form  of  a rapid  increase  in  pressure  at  about  the 
same  time  followed  by  an  abrupt  pressure  drop  at  a slightly  later  time, 
as  for  PS1  (Fig.  4.5a).  At  later  times  (after  about  the  3.3  msec)  the 
flow  pattern  is  more  like  a quasi-steady  flow  where  the  shock  tube  exit 
flow  can  be  considered  similar  to  a nozzle  or  jet  flow.  More  specifically, 
from  our  two-dimensional  REFLECT2  theoretical  calculations  (Appendix  A), 
it  appears  that  the  rapid  change  in  the  blast  pressure  (at  about  3.3  msec) 
corresponds  generally  to  the  arrival  of  the  contact  surface  between  the 
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"hot"  tunnel  gas  and  the  "cold"  gas  jet  from  the  shock  tube.  In  most 

figures  hereafter,  the  time  of  blast  arrival  at  the  transducer  and  a time 
* 

3.3  msec  later  are  indicated  as  a pair  of  vertical  arrows  in  order  to 
permit  the  reader  to  focus  his  attention  more  easily  on  the  range  of 
primary  significance  for  nuclear  blast  simulations  (between  the  arrows). 
While  the  test  results  for  larger  times  may  also  be  of  some  significance 
for  blast  simulations,  this  remains  to  be  demonstrated. 

Time  histories  of  ramp  and  cowl  pressures  are  presented  in  Figures 
4.5b  and  4.5c.  (See  Fig.  4.4  and  Table  4.1  for  transducer  locations). 

It  may  be  noted  that  there  is  considerable  variation  of  waveform  for  the 
different  locations  for  reasons  to  be  discussed  later.  (It  should  be 
noted  that  transducer  1990  generally  appears  to  be  truncated,  1904  is 
generally  unreliable  and  the  calibration  constant  used  for  2902  appears 
generally  unreliable) . 

Typical  time  histories  of  blast-induced  total  pressure  at  the 
blastward  engine  face  location  are  shown  in  Figure  4.5d.  Roughly,  the 
same  variations  are  seen  to  be  experienced  for  all  locations  and  the 
general  trends  are  similar  to  those  of  the  input  blast  pressure 
(e.g.,  PS2  in  Figure  4.5a). 

Typical  time  histories  of  blast  total  pressure  at  the  leeward  engine 
face  location  are  shown  in  Figure  4.5e.  The  leeward  pressures  are  somewhat 
similar  to  those  for  the  blastward  side  (Fig.  4.5d)  but  are  clearly 
different  from  those  for  the  blastward  side  in  two  significant  respects. 
First  the  initial  pressure  rise  is  more  rapid  for  the  leeward  side 
(essentially  a single  shock)  whereas  for  the  blastward  side  the  rise  time 
is  relatively  large  (say  about  0.3  msec),  and  generally  appears  to  involve 
a sequence  of  several  small  shocks.  Secondly,  the  blast  pressure  appears 
to  decay  more  rapidly  for  the  leeward  side  at  late  times. 

* 

For  a few  firings  the  onset  of  this  cold  gas  effect  appeared  to  occur 
somewhat  later,  but  for  most  cases  where  a definite  onset  could  be 
estimated,  it  appeared  to  occur  at  between  about  3.2  and  3.5  milliseconds 
after  blast  intercept. 
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Figure  4.5f  presents  times  histories  of  the  already  presented  three 
blast  pressure  measurements  PS1,  PS2  and  PS3  as  a fraction  of  pre-blast 
free  stream  total  pressure  (p^  ) and  also  presents  the  average  engine 
face  total  pressures  R20  and  R?I,  for  the  outboard  (blastward)  and  inboard 
(leeward)  inlets,  respectively.  As  would  be  expected,  these  last  two 
variables  correspond  closely  to  the  corresponding  variations  for  the 
individual  transducers  in  Figures  4.5d  and  4.5e,  respectively. 

Typical  claw  probe  results  are  shown  in  Figure  4.5g  and  4.5h. 

Directly  measured  parameters  are  the  "static"  pressure  PSi  and  the  left 
and  right  claw  arm  pressures,  PCLi  and  PCRi,  where  i designates  the 
probe  number.  Derived  quantities  are  the  resultant  flow  (sideslip) 
angle,  ANG,  the  local  Mach  number.  Mi,  and  the  total  pressure,  PTi. 

4.6  OBSERVATION  ON  BLAST  INPUT  CONDITIONS 

Before  proceeding  to  an  analysis  of  the  significance  of  the  test 
results  in  Section  5 it  is  important  to  first  indicate  briefly  the 
source  and  reliability  of  the  blast  input  overpressure  (Ap)  and  intercept 
angle  (40  data  presented  in  Table  4.4  and  to  indicate  what  differences 
in  blast  input  and  inlet  response  behavior  may  be  attributed  to  the 
different  locations  of  the  three  firing  shock  tubes. 

4.6.1  Blast  Input  to  the  Inlets 

The  strength  and  orientation  of  the  initial  blast  wave  incident  on 
the  inlets  and  its  subsequent  time  history  variations  must  be  known  in 
order  to  evaluate  the  significance  of  any  test  firing.  However,  insomuch 
as  the  blast  waves  for  these  tests  were  not  simple  constant  strength 
plane  waves  it  is  not  a simple  task  to  deduce  the  blast  strength  from 
the  test  data.  More  specifically,  the  orientation  of  the  blast  wave  at 
the  inlet  depends  on  the  tunnel  flow  conditions  (Mach  number),  on  the 
location  of  the  shock  tube,  and  on  the  shock  tube  firing  conditions 
(tube  driver  pressure/ tunnel  pressure)  and  may  be  time  dependent.  The 
intensity  of  the  blast  depends  on  the  radial  distance  and  polar  angle 
of  the  inlet  from  the  shock  tube  axis  and  on  the  tunnel  flow  conditions 
and  the  shock  tube  firing  conditions.  Furthermore  the  strength  of  the 
blast  wave  striking  the  inlet  and  the  blast  measurement  probes  can  be 
significantly  modified  by  interferences  between  the  blast  wave,  the 
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model  wing  and  fuselage  and  the  inlet,  particularly  at  late  times  after 
arrival  of  the  blast  at  the  inlet. 


4. 6. 1.1  Blast  Strength 

A first  approximation  for  the  blast  input  to  the  inlet  is  given  by 
the  static  pressure  transducers  mounted  on  the  claw  probes  as  indicated 
in  Figure  4.2.  Typical  time  histories  of  the  blast  signal  seen  by  these 
static  transducers  are  shown  in  Figure  4.6  for  firings  from  the  three 
shock  tubes. 

Referring  first  to  the  firing  from  tube  No.  1 (Fig.  4.6a),  probe  1 
(PS1)  is  struck  first  by  the  blast  wave,  probe  2 second,  and  probe  3 
last.  Probe  1 defines  the  blast  wave  as  it  is  seen  from  a position 
slightly  downstream  of  the  inlet  cowl  lip  and  about  6 inches  outboard 
of  the  cowl  lip.  This  record  should  provide  a reasonable  estimate  of 
the  incident  blast  strength  at  the  probe  location  for  a time  duration 
of  about  one  millisecond,  at  the  end  of  which  time  the  blast  wave  has 
reflected  from  the  blunt  side  of  the  inlet  model  and  returned  to  the 
probe,  thereby  providing  the  second  shock  pressure  increase  observed  on 
that  probe  at  that  time.  In  view  of  this  reflection  effect  it  is  difficult 
to  determine  the  late  time  history  of  the  incident  blast  wave  for  tube  1 
firings  from  probe  1 measurements.  (The  same  argument  applies  to  probe  3 
measurements  for  firings  from  tube  3,  as  can  be  seen  from  Figure  4.6c). 

The  blast  strength  time  history  obtained  from  probe  2,  located  above 
and  near  the  forward  end  of  the  inlet  appears  to  generally  give  the  best 
representation  of  the  blast  wave  input  that  can  be  obtained  without  a 
detailed  analysis  of  the  claw  probe  data.  There  are  no  obvious  strong 
interference  reflections  in  the  probe  2 data  (at  least  up  to  about  3.3 
milliseconds  after  blast  arrival),  although  there  certainly  must  be  some 
interferences  due  to  diffraction  of  the  blast  wave  about  the  inlet 
structure. 

The  blast  strength  record  obtained  from  probe  3 for  firings  from 
tube  1 or  2 (or  from  probe  1 for  tube  3 firings)  is  considered  only 
qualitatively  representative  of  the  true  blast  input  to  the  inlet  because 
the  blast  wave  experiences  substantial  diffraction  about  the  inlet 
structure  before  reaching  this  probe. 
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In  view  of  the  above  observations  and  in  view  of  frequent  malfunctions 
of  many  of  the  claw  probe  elements,  it  was  concluded  that  the  best  repre- 
sentations of  the  blast  strength  that  is  available  from  the  test  data  at 
present  are  the  probe  2 static  data  and  these  data  were  generally  used  to 
define  the  (nominal)  blast  strength  data  presented  in  the  tables  and 
figures  of  this  report.  A more  accurate  evaluation  of  the  true  incident 
blast  strength  and  its  time  history  derived  from  the  claw  probe  data 
would,  of  course,  be  desirable,  but  insufficient  time  was  available  to 
develop  this  subject. 

4. 6. 1.2  Blast  Orientation 

The  blast  incidence  angle,  <f>,  at  which  the  blast  wave  strikes  the 
inlet  is  of  considerable  importance  in  determining  blast  effects  on  the 
inlet.  The  blast  intercept  angle  is  defined  as  the  angle  between  a line 
normal  to  the  blast  front  and  the  axis  of  the  inlet,  with  0°,  90°  and  180° 
representing  head-on,  side-on  and  tail-on  blast  intercepts,  respectively. 

Initial  estimates  of  blast  intercept  angles  for  the  16T  tests  were 

made  from  the  shadowgram  studies  made  in  the  IT  tunnel  (Sec.  2.2).  Attempts 

were  made  to  take  similar  shadowgrams  during  the  16T  tests,  but  these 

were  generally  unsuccessful.  Instead,  estimates  of  the  blast  intercept 

angles  for  the  16T  test  conditions  were  generally  obtained  by  observing 

the  differences  in  arrival  times  of  the  blast  wave  at  claw  probe  2 and 

* 

at  the  transducers  1935,  1902  and  2902  located  near  the  inlet  mouth, 
under  the  assumption  that  the  blast  wave  could  be  considered  to  be 
locally  a plane  shock  wave,  travelling  with  respect  to  the  local  pre-blast 
fluid  at  a shock  speed  corresponding  to  the  shock  overpressure  in  Table 
4.4.  The  resulting  intercept  angle  estimates  are  presented  in  Table  4.4. 
Generally  these  angles  are  in  reasonable  agreement  with  estimates  from 
the  IT  test  shadowgrams,  for  cases  of  comparable  test  conditions. 


See  Figure  4.4. 


4.6.2  Shock  Tube  Location  Effects  on  Blast  Input  and  Inlet  Response 

The  blast  waves  produced  at  the  inlet  by  firings  from  different  shock 
tubes  can  produce  somewhat  different  effects  on  the  inlet  response  due  to 
factors  other  than  obvious  differences  in  intercept  angles  and  overpressures. 
It  is  the  purpose  of  this  section  to  clarify  this  question  somewhat  by 
examination  of  some  of  the  test  data  before  proceeding  to  a more  extensive 
analysis  of  the  data  in  Section  5. 

Consider  the  effects  of  firing  at  the  inlet  from  opposite  sides  of 
the  tunnel.  The  blast  waves  fired  at  the  model  from  shock  tubes  1 and  2 
encounter  the  inlet  essentially  without  any  initial  interferences  from 
the  model  wing  and  fuselage  body  (see  Fig.  3.1),  whereas  the  blast  wave 
from  shock  tube  3,  which  must  cross  the  wing  before  reaching  the  inlet, 
may  be  substantially  modified  before  striking  the  inlet.  Consequently, 
some  differences  might  be  expected  for  the  inlet  response  conditions  to 
firings  from  the  different  sides  of  the  tunnel.  In  order  to  indicate  the 
extent  of  these  differences,  Figure  4.7  presents  a comparison  of  inlet 
engine  face  total  pressure  time  histories  for  firings  from  shock  tubes 
2 and  3 for  roughly  comparible  overpressure  and  blast  intercept  conditions. 
The  left-hand  side  of  this  figure  presents  inlet  pressure  for  a firing 
from  tube  2 and  the  right  hand  side  presents  the  corresponding  inlet 
pressure  for  a firing  from  tube  3.  For  example,  transducer  1806  (in 
the  outboard  inlet)  for  a tube  2 firing  would  be  expected  to  have  the 
same  type  of  response  as  the  symmetrically  located  transducer  2831  (in 
the  inboard  inlet)  would  have  for  a firing  from  tube  3. 

It  is  evident  from  Figure  4.7a  that  the  inlet  response  pressures 
are  about  the  same  for  firings  from  the  two  tubes,  certainly  for  the 
estimated  duration  of  the  blast  event,  say  to  about  3.3  msec  after  blast 
arrival.  For  much  later  times,  there  are  appreciable  differences,  but 
this  is  not  of  great  importance  for  the  present  study. 

For  the  leeward  inlet,  Figure  4.7b  indicates  the  same  general  shape 
of  pressure  histories  for  the  two  firings  but  does  indicate  an  appreciably 
more  rapid  decay  of  pressure  for  the  firing  from  tube  3.  This  more  rapid 
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decay  can  be  attributed  partly  to  wing-fuselage- inlet  interference,  which 
would  be  expected  to  have  more  effect  on  the  leeward  inlet  than  on  the 
blastward  inlet,  and  partly  to  the  10°  difference  in  intercept  angles  for 
the  two  firings. 

Some  differences  are  also  observable  between  the  blast  inputs  and 
inlet  responses  for  firings  from  the  two  tubes  located  on  the  same  side 
of  the  tunnel.  For  example.  Figure  4.8  compares  inlet  ramp  pressure  time 
histories  for  firings  at  similar  overpressure  levels  from  the  side-by-side 
shock  tubes  1 and  2 (see  Fig.  4.4  for  transducer  locations).  It  may  be 
noted  that  the  two  sets  of  time  histories  are  quite  similar  during  the 
blast  event  (between  the  two  arrows) , especially  for  the  first  two  milli- 
seconds after  blast  arrival.  Then  the  pressures  tend  to  decrease  for  the 
rest  of  the  blast  event  for  the  tube  2 firing  and  to  increase  for  the  tube 
1 firing,  and  after  the  end  of  the  blast  event  the  two  sets  of  time 
histories  are  considerably  different.  The  greater  pressures  and  fluctua- 
tions observed  toward  the  end  of  and  after  the  blast  event  for  the  tube  1 
firing  may  be  attributed  partly  to  its  much  lower  blast  incidence  angle, 

84°  here  compared  to  110°  (see  Sec.  7),  and  partly  to  the  greater  closeness 
of  the  jet  from  tube  1 to  the  inlet  mouth. 


SECTION  V 

ANALYSIS  OF  16T  TEST  RESULTS 

This  section  presents  a general  discussion  and  analysis  of  some  of 
the  principal  effects  relevant  to  and  observed  from  the  16T  blast  test 
results.  First  a brief  discussion  is  given  of  the  basic  characteristics 
of  a blast  wave  interacting  with  an  inlet-engine  system  (Sec.  5.1), 
followed  by  a more  specific  discussion  of  the  transient  shock  wave 
pattern  in  an  inlet  after  blast  encounter  (Sec.  5.2).  Next  are  discussed 
distortion  effects  in  the  inlet  (Sec.  5.3).  Finally  some  observations 
are  made  on  the  ability  of  the  inlet  to  endure  long-duration  blast 
effects  (Sec.  5.4). 

5.1  GENERAL  PHENOMENA 

The  physical  characteristics  of  a blast  wave  and  its  initial  inter- 
actions with  an  inlet-engine  system  are  discussed  below.  The  initial 
discussion  deals  primarily  with  the  case  where  the  blast  wave  is  considered 
as  a simple  plane  shock  wave  of  constant  or  very  slowly  decaying  strength, 
resembling  the  blast  wave  from  a nuclear  weapon.  Subsequently,  attention 
is  given  to  the  differences  in  phenomena  which  are  associated  with  the 
more  complex  blast  wave  produced  by  the  shock  tubes  as  used  in  the 
present  studies. 

As  a blast  wave  travels  toward  an  aircraft  moving  in  space  it 
produces  changes  in  the  ambient  atmosphere  behind  its  front  which  are 
imposed  on  the  aircraft  as  it  is  enveloped  by  the  blast  field.  The 
blast  wave  is  defined  essentially  by  two  characteristics,  the  overpressure, 
Ap,  and  the  intercept  angle,  4> , between  a normal  to  the  blast  front  and 
the  direction  of  aircraft  motion.  With  reference  to  the  aircraft,  the 
flow  field  behind  the  blast  front  will  have  a greater  ambient  pressure, 

Po  + Ap,  a greater  total  pressure,  p^  + Ap^,  and  a blast-induced  increased 
side-slip  angle,  lj»  = + Alp.  Figures  5.1  and  5.2  illustrate  the 

magnitude  of  these  blast-induced  side-slip  angle  and  total  pressure 
changes  for  several  Mach  numbers.  For  example,  it  is  seen  from 
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(a)  MACH  0.55 

Figure  5.1.  Variation  of  blast-induced  side-slip  angle 
with  intercept  angle  and  overpressure. 
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Figure  5.1.  Continued 
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Figure  5.2.  Continued 
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Figure  5.2  Continued 
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Figure  5.1  that  the  greatest  blast-induced  yaw  angles  are  produced  by 
intercept  angles  somewhat  over  90°  (i.e.,  slightly  from  the  rear). 
Blast-induced  total  pressure  changes  in  Figure  5.2  are  seen  to  be 
largest  for  zero  intercept  angle  and  to  decrease  subsequentially  to 
relatively  small  values  for  tail-on  (<t>=180°)  intercepts. 

As  a blast  wave  envelopes  an  aircraft  these  blast  properties  may 
be  substantially  modified  by  diffraction  and  reflections  from  various 
parts  of  the  structure  before  and  after  the  wave  encounters  the  engine 
inlet . 

As  the  wave  proceeds  into  an  engine  inlet  it  may  experience  a 
sequence  of  reflections  from  the  walls  of  the  inlet  (as  discussed  in 
Sec.  5.2)  and  from  engine  components,  and  may  penetrate  in  part  completely 
through  the  engine  system. 

Eventually,  if  the  blast  wave  is  of  sufficient  duration,  these  blast 
reflection  effects  will  be  resolved  and  the  inlet-engine  system  will  come 
to  a quasi-steady  equilibrium  condition  consistent  with  the  blast-induced 
overpressure,  total  pressure  and  side-slip  angle,  characterized  by  an 
engine  mass  flow  rate  which  can  be  different  from  that  for  the  pre-blast 
condition . 

In  discussing  the  effects  of  a blast  wave  on  inlet-engine  performances 
it  is  convenient  to  think  of  the  phenomena  in  terms  of  the  following 
somewhat  overlapping  three  stages  of  the  blast  encounter. 

The  first  or  initial  stage  consist  of  the  time  during  which  the 
blast  wave  moves  down  the  inlet  into  the  engine,  reflects  in  part  from 
the  engine  structure  or  fan  of  a turbofan  engine  and  subsequently  moves 
upstream  and  out  the  mouth  of  the  inlet.  In  the  case  of  the  present 
0.1-scale  model  the  time  required  for  these  events  up  to  the  time  when 
the  reflected  shock  wave  has  moved  upstream  of  the  simulated  engine 
face  section  is  on  the  order  of  2 milliseconds.  (Time  zero  is  defined 
here  as  the  time  when  the  blast  wave  first  strikes  the  cowl  lip  of  the 
blastward  inlet.) 

The  second  stage  of  the  blast  flow  encounter  is  associated  with  the 
inception  and  convection  of  boundary-layer  disturbances  down  the  inlet. 

In  particular,  boundary- layer  separation  and  vortex  formations  may  be 
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produced  at  the  cowl  lip  by  large  blast-induced  side-slip  angles  and  will 
be  convected  downstream  at  speeds  comparable  to  or  lower  than  the  particle 
velocity  in  the  inlet,  which  is  generally  much  smaller  than  the  rate  at 
which  a blast  wave  moves  down  the  inlet.  In  the  case  of  the  tested  inlet, 
it  would  take  on  the  order  of  3 milliseconds  (after  cowl  lip  blast 
encounter)  for  such  particle  velocity  effects  to  reach  the  engine  face. 
Other  boundary  layer  effects  which  might  be  expected  in  this  same  time 
period  would  be  boundary  layer  instabilities  excited  by  the  reverse  flow 
reflected  upstream  from  the  engine  throat.  For  the  present  inlet,  such 
effects  would  not  be  expected  to  begin  at  the  engine  face  before  times  of 
about  2 milliseconds  after  blast  intercept. 

The  third  stage  of  the  blast  flow  consists  of  the  late  time  stage 
when  all  significant  internal  shock  wave  systems  have  been  dissipated  or 
reflected  from  the  inlet-engine  system  and  boundary  layer  effects  have 
reached  a quasi-steady  equilibrium  condition.  This  stage  was  not  reached 
in  the  present  tests  because  of  the  limited  duration  of  the  blast  waves 
used . 

5.1.1  Differences  Between  Shock  Tube  and  Nuclear  Blast  Flows 

The  initial  flow  from  the  shock  tubes  used  in  the  present  tests 
generally  resembles  the  essential  features  of  a nuclear  blast  wave  at 
least  for  part  of  the  blast-type  flow  period  described  (in  Section  4.5)  of 
about  3.3  milliseconds.  However  there  are  some  noteworthy  features, 
which  should  be  borne  in  mind,  where  the  shock  tube  flow  differs  from 
that  for  a nuclear  blast.  First  the  strength  of  the  shock  tube  blast 
wave  remains  relatively  constant  for  only  one  or  two  milliseconds,  after 
which  its  strength  generally  tends  to  decrease  toward  zero  at  4 to  6 
seconds  after  initial  shock  arrival.  Also,  in  the  late  decay  stage, 

"cold"  air  blown  out  of  the  shock  tube  produces  a high  velocity  jet  with 
a complex  spatial  variation  of  fluid  properties  which  bears  little 
resemblance  to  a nuclear  blast  field.  As  a general  rule  of  thumb  for 
the  present  tests,  it  appears  that  this  "cold"  jet  does  not  strongly 
influence  the  test  results  for  times  after  blast  arrival  less  than  about 
3.3  milliseconds  (see  Sec.  4.5). 
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A sketch  illustrating  the  above-discussed  sequence  of  events  in  the 
blastward  inlet  is  presented  as  Figure  5.3  in  the  form  of  a time  iistance 
plot  for  an  inlet  mass  flow  rate  of  350  lb/sec.  The  time  origin  is  the 
time  when  the  blast  wave  first  arrives  at  the  cowl  lip.  The  initial 
blast  wave  moves  down  the  inlet  to  the  simulated  engine  throat  vane 
system  and  is  reflected  back  upsteam  as  indicated  by  the  "shock"  curve. 
The  "particle  velocity"  curve  indicates  the  convective  speed  of  material 
in  the  inlet;  "vortex"  motions  are  indicated  here  as  travelling  at 
speeds  between  1/2  and  1 times  the  particle  velocity  and  are  indicated 
to  reach  the  engine  face  at  times  on  the  order  of  3-6  milliseconds, 
depending  on  the  vortex  speed.  The  cold  gas  motion  is  indicated  by  the 
two  upper  curves,  starting  at  time  3.3  msec.  This  cold  gas  effect  might 
be  expected  to  travel  down  the  inlet  at  some  speed  between  the  particle 
velocity  and  the  sound (+  particle)  velocity,  depending  on  how  close  the 
inlet  mouth  is  to  the  jet  from  the  shock  tube.  It  is  important  to  note 
here  that  large  shock-type  cold  gas  effects  (see  lower  cold  gas  curve) 
can  arrive  at  the  engine  face  before  vortex  effects  arrive,  hence  making 
it  difficult  to  unambiguously  evaluate  vortex  effects  at  the  engine  face 
by  simple  inspection  of  test  records. 

5.2  SHOCK  PATTERN  IN  THE  INLETS 

As  a guide  to  the  discussion  in  later  sections  of  blast-induced 
distortion  and  its  effects  on  engine  performance,  this  section  briefly 
discusses  some  details  of  the  character  of  the  blast-induced  shock  flow 
pattern  in  the  inlet. 

Figure  5.4  presents  a sketch  of  the  tested  inlet  configuration  and 
shows  the  shock  wave  pattern  inside  the  inlet  for  a time  when  the  blast 
shock  has  penetrated  down  the  blastward  (outboard)  inlet  to  slightly 
beyond  the  engine  face  location.  In  this  figure,  the  blast  is  assumed 
to  have  struck  the  inlet  about  perpendicular  to  the  inlet  axis  (<J>=90°) 
from  the  outboard  side  of  the  inlet.  After  encountering  the  lip  of  the 
outboard  cowl,  the  blast  shock  wave  begins  to  diffract  about  the  cowl 
lip  and  to  penetrate  into  the  inlet  at  a downstream  speed  about  equal 
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figure  5.3.  Inlet  time-distance  relationships . 
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Figure  5.  A.  Typical  shock  wave  pattern  in  the  inlets. 
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to  the  sum  of  the  (pre-blast)  flow  velocity  in  the  inlet  plus  the  speed 
of  sound  in  the  inlet.  At  the  same  time  the  part  of  the  shock  wave 
upstream  of  the  cowl  lip  proceeds  to  the  outboard  ramp  surface,  experiences 
essentially  normal  reflection  there,  with  about  doubling  of  its  intensity, 
and  then  tends  to  bounce  back  and  forth  between  the  ramp  and  cowl  surfaces 
while  being  at  the  same  time  convected  downstream  into  the  inlet.  The 
resulting  overall  blast  front  pattern  at  any  time  can  be  constructed 
graphically  by  assuming  the  shocks  to  travel  at  essentially  the  speed  of 
sound  relative  to  the  local  ambient  inlet  velocity  and  the  accuracy  of 
the  construction  can  be  checked  or  improved  by  adjusting  these  constructions 
slightly  to  conform  to  the  experimentally  observed  times  of  blast  arrivals 
at  the  various  transducer  stations  along  the  ramp  and  cowl  indicated  in 
Figure  5.4.  (Generally  the  first  two  or  three  shock  arrival  times  can 
be  clearly  distinguished  in  the  transducer  time  histories.)  It  should 
be  noted  that  the  shock  pattern  in  Figure  5.4  is  rather  simplified  in 
that  it  assumes  regular  reflection  of  the  waves  at  the  cowl  and  ramp 
surfaces  and  it  does  not  consider  the  interference  effects  of  the  engine 
hub  structure  (Figure  4.1).  Actually  Mach  stems  will  tend  to  form  for 
some  reflected  waves  and  additional  reflections  will  occur  from  the  hub; 
however.  Figure  5.4  does  illustrate  the  overall  pattern  well  enough  for 
the  present  discussions. 

The  shock  pattern  on  the  leeward  side  of  the  inlet  is  much  simpler. 
Here,  the  shock  wave  simply  diffracts  about  the  tip  of  the  ramp  and 
tends  to  enter  the  leeward  inlet  essentially  as  a one-dimensional 
disturbance  having  a single  shock  with  its  front  essentially  normal 
to  the  inlet  side  walls. 

The  shock  patterns  shown  in  Figure  5.4  were  based  partly  on  test 
data  for  Run  39  (Part  544)  conditions,  Mach  0.85  at  350  lb/sec,  but 
would  be  expected  to  apply  rather  well  to  any  of  the  other  test  Mach 
numbers  as  well. 

Now,  to  indicate  the  significance  of  the  type  of  shock  pattern  seen 
in  Figure  5.4,  consider  first  the  blastward  inlet.  It  is  seen  that  the 
shock  waves  are  generally  not  at  all  parallel  to  the  inlet  walls,  except 
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for  the  initial  blast  front.  Hence,  there  is  generally  an  appreciable 
time-varying  spatial  variation  of  blast  pressure  across  any  section  of  the 
inlet,  which  lasts  from  the  time  of  first  blast  arrival  until  that  time 
when  all  of  the  incident  blast  front  and  its  reflections  will  have  been 
either  convected  downstream  through  the  engine  or  reflected  from  the 
engine  throat  back  upstream  and  out  the  mouth  of  the  inlet.  This  spatial 
pressure  variation  across  any  section  of  the  inlet  will,  of  course, 
produce  flow  distortion  in  the  section,  the  intensity  of  which  will  be 
related  to  the  intensity  of  the  shock  or  shocks  passing  through  or  near 
to  that  section  at  the  time  in  question. 

Another  consequence  of  the  complex  shock  pattern  in  the  blastward 
inlet  is  its  effect  on  the  rate  at  which  the  shock  pressure  builds  up. 
Essentially,  for  the  blastward  inlet,  the  pressure  buildup  consists  of  a 
sequence  of  separated  shocks  adding  to  each  other,  so  that  the  resultant 
effect  is  a relatively  slow  staircase  type  buildup  of  pressure,  taking 
on  the  order  of  0.3  millisecond  for  typical  cases.  On  the  leeward 
inlet,  however,  where  there  is  essentially  a single  shock  front,  the 
pressure  buildup  is  essentially  instantaneous. 

5.3  ANALYSIS  OF  INLET  DISTORTION 

A primary  objective  of  the  test  program  was  to  obtain  time  histories 
of  the  inlet  distortion  parameters  resulting  during  blast  wave  inlet 
interactions.  In  this  study,  time  histories  of  the  following  inlet 
distortion  parameters  were  calculated  from  the  test  data  for  both  inlets 
for  all  runs  and  were  plotted  as  illustrated  in  Figure  5.5  for  Run  8 
(Part  573):  IDCi,  IDRi,  IDC12,  IDC45,  IDR,  IDC,  IDL,  IDA  and  IDT,  where 


i designates  the  ith  instrumented  engine  face  ring  (i  = 1 to  5) . These 
particular  distortion  parameters,  defined  in  Table  5.1,  were  selected 
because  they  are  standard  for  the  analysis  of  the  B-l  inlet  and  engine 
system.  The  discussion  of  distortion  in  this  report  is,  however, 
restricted  to  the  parameters  IDC,  IDR  and  IDL,  and  primarily  to  IDL.  IDC 
is  the  total  circumferential  distortion  index,  IDR  is  the  total  radial 
distortion  index  and  IDL  is  the  fan  stall  margin  ratio.  IDL  is  an 
overall  measure  of  the  inlet  distortion,  which  depends  on  IDC, 
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Figure  5.5.  Concluded. 


16.00 


20.00 


TABLE  5.1 


equations  for  calculating  distortion 

PARAMETERS  IDC,  IDR,  IDL,  IDA  AND  IDT 


Circumferential  Distortion  on  Individual  Rings 


IDCi  " <Pt2 


“ P 


avci 


t2 


)/p 


min . 


t2 


where 


i = 1 to  5 rings 

Hub  Distortion 

IDC12  = (IDC1  + IDC2)/2 

Tip  Distortion 


IDC45  = (IDC4  + IDC5)/2 
Total  Circumferential  Distortion 


IDC  = Largest  of  IDC10  or  IDC. _ 

Radial  Distortion  on  Individual  Rings 

Pt2 

IDRi  = Largest  of  (1.0  - avej  ,0) 

Pt2 

Total  Radial  Distortion 

IDR  = Largest  of  IDR.^  IDR2>  IDR4>  or  IDR5 
Fan  Stall  Margin  Ratio 

IDL  = b (KCIRC) IDC  + (KRAD)IDR 

where 


b = 


IDR/ IDC 

A+B (IDR/ IDC)  + C 
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MACH 

NUMBER 

TUNNEL  TOTAL 
TEMPERATURE  (°R) 

1 

KCIR 

KRAD 

0 < .15 

All 

7.93 

12.41 

.15  < .35 

All 

7.69 

11.76 

.35  < .65 

All 

7.69 

12.28 

.65  < .75 

All 

7.69 

12.42 

.75  < 1.5 

Less  Than/Equal 

To  560. 

7.69 

12.83 

.75  < 1.5 

Greater  Than  560. 

7.93 

12.99 

1.5  < 2.3 

All 

7.69 

11.75 

Average  Distortion: 

IDA  = (p  - pt  )p 


t2  t . t2 
min 


Total  Distortion: 


IDT  = (pt  - pt  )/pt 
max  min 


SUPER  POSITION  FACTOR  b 


IDR  and  the  speed  altitude  regime.  An  indication  of  the  normally 
allowable  range  of  these  distortion  parameters  for  subsonic  flight  of 
the  B-l  aircraft  is  given  in  Figure  5.6  (from  Ref.  5.1).  Generally,  the 
inlet  and  engine  can  be  considered  to  perform  satisfactorily  at  values 
of  IDL  up  to  unity,  which  corresponds  to  values  of  IDC  and  IDR  to  the 
left  of  and  below  the  solid  line  curve  in  Figure  5.6.  For  larger  values 
of  IDC,  IDR  or  IDL  the  engine  can  stall,  but  will  not  necessarily  do  so. 

5.3.1  Data  Presentation 

Typical  time  histories  of  distortion  parameters  for  various  over- 
pressure levels  are  presented  in  Figures  5.7  through  5.17.  Figure  5.7 
presents  IDL  for  the  lower  tested  mass  flow  rate  (~300  lb/sec)  at  Mach 
0.70  for  both  inlets  and  Figures  5.8  and  5.9  present  similar  data  for 
the  higher  tested  mass  flow  rate  («350  lb/sec)  for  the  blastward  and 
leeward  inlets,  respectively.  Figure  5.10  presents  IDC,  IDR  and  IDL  for 
the  lower  mass  flow  rate  («300  lb/sec)  at  Mach  0.85  for  the  blastward 
inlet  and  Figure  5.11  presents  IDL  for  the  leeward  inlet;  Figure  5.12 
presents  IDC,  IDR  and  IDL  for  the  higher  mass  flow  rate  («350  lb/sec). 
Mach  number  effects  on  IDC,  IDR  and  IDL  are  indicated  in  Figures  5.13 
and  5.14  for  blastward  and  leeward  inlets  for  firings  from  Tube  2 and, 
for  IDL  only,  in  Figures  5.15  and  5.16  for  firings  from  Tube  1.  IDL 
measurements  made  during  the  three  yawed  firings  are  presented  in 
Figure  5.17. 

5.3.2  Effects  at  Low  Mass  Flow  Rates 

Considering  now  the  distortion  effects  at  the  lower  tested  mass 
flow  rate  of  about  300  lb/sec  (full  scale),  it  can  be  seen  from  Figure 
5.7  for  Mach  0.70  and  from  Figure  5.10c  and  5.11  for  Mach  0.85  that  the 
stall  margin,  IDL,  is  usually  well  below  the  normally  safe  limit  of 
unity,  usually  being  near  to  or  below  0.5,  for  the  time  range  where  the 
test  data  are  definitely  representative  of  a blast  wave  (between  the  two 
arrows  in  the  figures).  Corresponding  values  of  IDC  and  IDR,  shown  in 
Figure  5.10  for  Mach  0.85  are  similarly  rather  low,  being  generally 
below  0.10  for  IDC  and  below  0.05  for  IDR. 
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Figure  5.6.  Range  of  distortion  indices  for  B-l  inlet  at  subsonic  cruise 
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Figure  5.9.  Distortion  time  histories  for  three  overpressures  at  Mach  0.70 
for  the  leeward  inlet,  flow  rate%  350  lb/sec,  tube  2. 
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Figure  5.10.  Distortion  time  histories  for  four  overpressures  at  Mach  0.85 
for  the  blastward  inlet,  flow  rate  300  lb/sec,  tube  1. 


OUTBOARD  INLET 


C_) 

Q 


INBOARD  INLET 


(_> 

Q 


0.25 
0-20 
0- IS 
o.io 
0.05 
0.00 


i 1 i ( 1 1 

0-00  4.00  8-00  12-00  16-00  20-00 

TIME  IN  MILLISECONDS 

Figure  5.12.  Distortion  time  Histories  for  two  overpressures  at  Mach  0.85 
flow  rate w 150  lh/sec.  tube  2. 
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Figure  5.14.  Distortion  time  histories  for  four  Mach  numbers  for  the 
leeward  inlet,  flow  rate  <s  350  lb/sec,  tube  2. 


In  a few  cases,  values  of  IDL,  IDC  and  IDR  in  excess  of  normally 
allowable  limits  were  observed.  For  example,  for  Run  28  (Fart  584)  in 
Figure  5.10,  large  distortion  values  were  experienced  in  the  blastward 
inlet  shortly  after  blast  intercept,  but  the  duration  of  this  large 
distortion  appears  to  have  been  less  than  0.1  msec  (1  msec  full-scale). 

It  is  uncertain  whether  the  duration  of  this  distortion  pulse  was  long 
enough  to  imply  a significant  adverse  effect  on  engine  performance.  As 
a matter  of  possible  interest,  Figure  5.18  presents  engine-face  total- 
pressure  contours  for  this  run,  indicating  the  total-pressure  variation 
across  the  engine  face  at  a time  close  to  the  time  of  peak  distortion 
(indicated  by  the  arrow  in  the  inset  IDL  plot) . The  corresponding 
instantaneous  values  of  IDL,  IDC  and  IDR  are  1.05,  0.103  and  0.061, 
respectively.  The  numerical  values  in  this  figure  represent  the  difference 
between  the  face-average  total-pressure  and  the  local  total-pressure, 
expressed  as  a fraction  of  face-average  total-pressure.  Test  data  points 
(or  interpolations  thereof  for  a few  bad  data  points)  are  represented  as 
large  numbers  and  contour  levels  by  small  handwritten  numbers.  It  may 
be  noted  that  the  lowest  total-pressure  drop  is  about  17  percent  below 
the  face-average  total-pressure. 

5.3.3  Effects  at  High  Mass  Flow  Rates 

The  distortion  effects  at  the  higher  tested  mass  flow  rate  of  350 
lb/sec  are  generally  at  least  somewhat  larger  than  for  the  lower  mass 
flow  rate.  However,  for  the  lower  overpressures,  up  to  about  4 psi,  the 
stall  margin,  IDL,  is  generally  observed  to  be  well  below  the  limit  of 
unity  in  the  blast  range  (between  the  two  arrows),  being  generally  near 
to  or  below  0.9  (see  Figures  5.8  and  5.9  for  Mach  0.70,  Figure  5.12  for 
Mach  0.85  and  Figures  5.13  and  5.14  for  a range  of  Mach  numbers).  For 
higher  overpressures  than  4 psi,  IDL  sometimes  appears  to  exceed  unity 
in  the  blast  range  (e.g.,  for  the  blastward  inlet  at  Ap  = 5.2  psi  at 
Mach  0.70  in  Figure  5.8  and  for  the  leeward  inlet  for  Ap=5.8  psi  at  Mach 
0.85  in  Figure  5.12),  but  only  for  a small  fraction  of  a millisecond  at 
most.  Outside  the  blast  range  (to  the  right  of  the  second  arrows  in 
the  figures)  IDL  values  appreciably  larger  than  unity  are  observed 
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for  some  firings  (e.g.,  for  the  blastward  inlet  at  Ap  = 3.2  psi  at  Mach 
0.7  in  Figure  5.8  and  for  blastward  and  leeward  inlets  at  Ap  = 3.8  psi 
at  Mach  0.55  in  Figure  5.13),  but  further  study  would  be  required  to 
verify  whether  these  late  time  data  are  relevant  to  the  blast-inlet 
problem.  Possible  sources  of  late-time  distortion  are  discussed  in 
Section  9. 

5.3.4  Mach  Number  Effects 

No  definite  significant  effects  of  Mach  number  on  the  distortion 
indices  IDL,  IDC  and  IDR  were  observed  in  the  blast  range  (e.g. , see 
Figures  5.13  and  5.14).  This  lack  of  observed  effect  is  attributed 
primarily  to  the  choking  action  of  the  vane  system  downstream  of  the 
simulated  engine  face,  which  tends  to  maintain  the  internal  inlet  flow 
at  any  location  in  the  inlet  at  a constant  Mach  number  condition,  regard- 
less of  the  state  of  the  external  flow  conditions.  Further  more  definite 
evidence  indicating  this  apparent  lack  of  Mach  number  effects  on  inlet 
pressures  is  given  in  Figure  5.19,  which  compares  inlet  ramp  pressures 
for  transducer  1980  (located  about  half-way  down  the  inlet,  see  Figure 
5.4)  for  four  Mach  number  conditions.  It  may  be  noted  that  most  fine 
details  of  the  pressure  histories  in  Figure  5.19  are  almost  identical 
for  the  different  Mach  numbers,  especially  for  Mach  0.55  and  0.85. 

5.3.5  Effects  of  Yaw  Angle 

No  definite  effects  of  yaw  angle  on  distortion  are  readily  evident 
from  the  results  of  the  three  yawed  test  firings  (Figure  5.17).  Appreciable 
increases  in  IDL  occurred  for  both  inlets  in  each  firing,  but  the 
increases  did  not  generally  appear  much  greater  than  for  similar  unyawed 
firings  (see  Figs.  5.8  and  5.9).  However,  it  should  be  noted  that  a 
significantly  high  level  of  distortion  appears  to  have  developed  in  the 
leeward  inlet  toward  the  end  of  the  nominal  blast  event  (at  second  arrow) 
of  Run  25  (Fig.  5.17b). 

5.3.6  Long  Duration  Distortion  Effects 

If  the  rate  of  decay  of  a blast  wave  is  extremely  small,  as  for  a 
large  nuclear  burst,  then  the  struck  inlet-engine  system  will  tend  at 
late  times  to  come  to  a quasi-steady  equilibrium  condition  where  the 
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effect  of  the  blast  wave  will  be  essentially  the  same  as  would  be 
produced  by  steady  yawed  flight  at  the  side-slip  angle  and  total  pressure 
caused  by  the  blast  (as  given  by  Figures  5.1  and  5.2.  In  regard  to  this 
point,  it  may  be  noted  that  our  test  results  for  steady  yawed  conditions 
generally  indicate  10%  or  less  increase,  or  decrease,  in  the  distortion 
parameter  IDL  as  the  sideslip  angle  varies  from  the  zero  angle  condition 
to  angles  anywhere  in  the  range  from  -10°  to  +10°.  Hence,  it  might  be 
expected  that  late-time  blast  distortions  would  be  small  for  blast 
overpressures  producing  blast-induced  sideslip  angles  up  to  at  least  10°, 
as  given  by  Figure  5.1  in  Section  5.  For  example,  for  Mach  0.85  flight. 
Figure  5.1  indicates  that  blast  overpressures  up  to  about  3.3  psi  would 
not  produce  sideslip  angles  over  10°  and,  hence,  would  not  be  expected 
to  produce  large  late-time  distortion  effects. 

The  value  of  the  above  observations  is,  of  course,  of  limited  value 
since  we  cannot  at  this  time  specify  what  we  mean  by  "late  time"  in  a 
very  definite  manner.  However,  judging  from  BID  code  calculations  for 
times  up  to  about  6 milliseconds,  presented  later  in  Section  7 (Figure 
7.2),  it  appears  that  quasi-steady  equilibrium  would  not  be  reached  for 
the  test  model  until  times  considerably  greater  than  6 milliseconds  (or 
60  milliseconds  for  a full  scale  B-l  inlet) . 

5.3.7  Concluding  Remarks  on  Distortion  Effects 

In  concluding  this  discussion  of  distortion  effects,  it  should  be 
noted  that  while  generally  very  few  large  distortion  values  were  obtained 
during  the  early-time  definitely  blast-type  flow  periods  of  the  present 
tests,  the  distortion  criteria  used  above  (IDC,  IDR  and  IDL)  are  not 
necessarily  completely  reliable  indices  of  whether  an  engine  would  stall 
under  transient  blast  conditions.  Other  distortion  criteria  should  also 
be  considered  (e.g.,  Ref.  5.2  and  5.3)  and  attention  shoudl  also  be 
directed  to  the  actual  contours  of  pressure  distortion  which  are  obtained 
under  blast  intercept  conditions,  such  as  is  illustrated  in  Figure  5.18 
for  one  run  at  a time  shortly  after  blast  intercept  (see  arrow  in  inset 
IDL  plot)  when  appreciable  distortion  was  experienced  in  the  outboard 
inlet.  Also  attention  should  be  directed  to  the  durations  of  any  large 
observed  distortion  levels,  comparing  them  with  duration  results 


from  other  tests  (e.g..  Ref.  5.4).  Furthermore,  attention  should  be 
directed  also  to  effects  of  blast-induced  longitudinal  pressure  variations 
(Sec.  9)  and  blast-fan  interactions  (Sec.  8)  to  assess  more  completely 
the  possible  stall  inducing  potential  of  blast  waves. 

These  same  types  of  questions  should  also  be  considered  with  respect 
to  the  large  distortion  values  sometimes  observed  at  times  after  the 
fairly  definite  blast  type  flow  duration  time  of  3.3  milliseconds. 

Most  of  the  above  topics  are  not  explored  in  the  present  report 
because  the  reduced  test  data  were  not  made  available  in  a suitable  form 
in  time  to  permit  such  studies  to  be  performed. 
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SECTION  VI 

COMPARISON  OF  THEORY  AND  EXPERIMENT 


Concurrent  with  the  performance  of  this  experimental  investigation, 
Kaman  AviDyne  developed  a theoretical  two-dimensional  computer  code, 
designated  BID  (blast-induced-distortion) , for  predicting  the  transient 
flow  field  produced  in  a inlet  by  a blast  wave  striking  the  inlet  at 
an  arbitrary  angle  of  incidence  (Reference  6.1).  Calculations  of  inlet 
pressure  time  histories  were  made  with  this  code  for  four  conditions 
similar  to  those  of  the  16T  experimental  tests  and  the  results  of  these 
calculations  are  compared  with  the  test  results  in  this  section. 

The  specific  conditions  for  which  BID  calculations  were  made  are: 
constant  shock  overpressure  of  5 psi  at  one  atmosphere  ambient  pressure, 
inlet  mass  flow  rate  (full-scale)  350  lb/sec,  tunnel  Mach  number  of  0.70 
for  a blast  intercept  angle  of  <J)=90°  and  tunnel  Mach  number  of  0.85  for 
intercept  angles  of  4>=90°,  105°  and  135°.  Calculated  BID  pressure  time 
histories  for  various  locations  in  the  inlets  are  presented  in  Figures 
6.1  to  6.5  together  with  experimental  data  for  those  test  runs  corres- 
ponding most  closely  to  the  conditions  of  the  theoretical  calculations. 

The  reader  is  reminded  here  that  in  looking  at  the  experimental  data 
in  Figures  6.1  to  6.3,  attention  should  be  directed  on  times  between  the 
two  arrows  above  each  curve;  the  data  at  later  times  are  not  considered 
necessarily  representative  of  a blast  type  flow  (See  Sec.  4).  Also  it 
should  be  recognized  that  the  theoretical  calculations  shown  in  Figures 
6.1  to  6.4  are  for  a constant-strength  (step)  blast  wave  (5  psi  over- 
pressure) , whereas  the  test  data  correspond  to  a time-varying  blast  wave 
characterized  by  an  initial  shock  value  (shown  below  each  figure),  and 
possible  subsequent  increases  in  pressure,  followed  by  a tendency  to 
decrease  to  zero  in  4-6  milliseconds  after  initial  shock  arrival. 

The  reader  is  also  advised  to  make  comparisons  of  theory  and 
experiment  primarily  on  the  basis  of  wave  shape  rather  than  on  amplitude 
comparisons,  since  the  calibration  factors  used  for  some  of  the  shown 
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ramp  and  cowl  data  are  either  nominal  or  questionable  values  and  a few 

* 

appear  to  be  significantly  in  error  . 

Considering  first  the  Mach  0.7  condition,  it  can  be  seen  from  the 
presentation  of  ramp  and  cowl  time  histories  in  Figure  6.1  that  there  is 
a great  difference  in  the  pressure  waveforms  for  different  locations  in 
the  inlet  (see  Figure  5.4  to  identify  the  locations  indicated  in  the 
figure) , and  that  in  almost  every  case  the  theoretical  results  follow 
well  most  of  the  details  observed  in  the  experimental  data.  For  example, 
for  location  1970,  the  BID  calculation  clearly  delineates  the  first 
three  shock  arrivals  and,  at  late  time,  delineates  a large  pressure 
increase  resulting  from  reflection  of  the  blast  from  the  simulated 
engine  throat  of  the  inlet  model. 

A similar  comparison  is  made  for  ramp  pressures  for  Mach  0.85  in 
Figures  6.2  and  6.3  and  for  cowl  pressures  in  Figure  6.4.  Again  the 
theoretical  and  experimental  wave  shapes  are  generally  in  good  agreement. 

A few  noticible  differences  do  appear  to  exist,  such  as  the  more  rapid 
late  time  decays  (say  near  the  second  arrow)  of  the  experimental  pressures 
for  locations  1950  and  1970  in  Figure  6.3  and  for  the  cowl  locations  in 
Figure  6.4.  These  differences  may  be  attributed  in  part  to  differences 
between  the  constant  strength  blast  assumed  for  the  calculations  and  the 
actual  variable  strength  blast  of  the  test. 


A final  comparison  is  made  in  Figure  6.5  between  calculated  and 
experimental  total  pressures  at  the  simulated  engine  face  location  for 
a location  near  the  ramp  in  the  blastward  (outboard)  inlet.  The  two 
transducers  shown  are  located  vertically  above  one  another  in  the  inlet 
(see  Figure  4.3)  so  that  the  theoretical  (two-dimensional)  BID  prediction 
is  the  same  for  both  transducers.  Differences  between  the  two  experimental 
results  are  hence  indicative  of  the  importance  of  three-dimensional 
effects  not  covered  by  the  theory  in  its  present  form. 


These  calibration  factor  problems  arose  as  a result  of  pre-test  mechanical 
and  electrical  malfunctions  of  the  transducers  which  could  not  be  resolved 
in  the  very  limited  time  available  for  the  16T  tunnel  tests.  Preliminary 
examinations  of  the  test  data  indicated  that  more  reliable  estimates 
for  some  of  these  factors  can  be  obtained  by  considering  some  redundant 
features  of  the  data,  but  this  has  not  yet  been  done. 
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Comparison  of  theoretical  and  experimental  time  histories  of  ramp 
pressures  in  the  blastward  inlet  at  Mach  0.85,  tube  1 data, 
flow  rate  «=  <50  lb/ sec. 
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Figure  6.3.  Comparison  of  theoretical  and  experimental  time  histories  of  ramp 
pressures  in  the  blastward  inlet  at  Mach  0.85,  tube  2 data, 
flow  rate  « 350  lb/sec. 


In  order  to  permit  a reasonable  quantitative  comparison  of  theory 
and  experiment  in  Figure  6.5,  it  was  necessary  to  take  into  account  the 
fact  that  the  pre-blast  (steady-state)  pressures  were  slightly  different 
and  the  input  blast  overpressures  and  orientations  were  significantly 
different  for  the  two  cases.  This  difference  was  taken  into  account  by 
multiplying  the  theoretical  BID  code  curve  by  the  ratio  of  the  experimental 
and  theoretical  values  of  the  shock  jump  in  input  total  pressure  and  by 
vertically  shifting  the  theoretical  curve  so  that  the  pre-blast  steady- 
state  values  were  the  same  for  theory  and  experiment. 

The  theoretical  and  experimental  total  pressures  in  Figure  6.5  are 
seen  to  be  in  reasonably  good  agreement  in  the  following  respects.  Both 
indicate  a reasonably  similar  (slow)  rate  of  rise  of  the  pressure  to  a 
peak  value  in  about  0. 3-0.4  millisecond  and  both  have  about  the  same 
mean  peak  level  of  pressure  for  about  2 milliseconds.  Also  the  third 
peak  in  the  theoretical  curve  (at  about  1 msec  after  blast  arrival), 
which  represents  the  return  of  the  blast  wave  reflected  from  the  simulated 
engine  throat,  is  in  agreement  with  the  corresponding  experimental  time. 
There  are,  of  course,  some  noticable  differences,  in  that  the  experimental 
peak  values  are  appreciably  greater  than  the  theoretical  values  and,  at 
later  times,  may  be  smaller.  These  differences  may  be  attributed  in 
part  to  the  finite  resolution  of  the  numerical  code,  some  attenuation  of 
the  shock  for  the  code  due  to  an  effective  "artificial  viscosity",  with 
a fixed-mesh  cell  system,  to  three-dimensional  effects  (as  evidenced  by 
differences  between  the  two  experimental  traces)  not  covered  by  the 
theory,  and  to  the  use  in  the  theoretical  calculation  of  the  approximation 
that  the  blast  wave  is  of  constant  strength. 

In  summary,  the  above  limited  comparisons  indicate  that  the  major 
features  of  the  transient  pressures  observed  in  the  inlet  tests  are  well 
represented  by  the  BID  code  results.  However,  a more  extensive  compara- 
tive study  of  the  16T  results  would  be  desirable  to  more  definitively 
establish  the  limitations  of  the  code  predictions.  Such  a correlation 
would  require  incorporation  into  the  BID  code  of  the  actual  time  history 
of  the  blast  wave  incident  on  the  inlet  (instead  of  the  constant  over- 


130 


pressure  value  used  here).  This,  in  turn,  requires  evaluation  of 
the  input  blast  time  history  from  the  test  data,  which  is  a non- 
trivial task  because  of  possible  significant  fuselage  inlet  and 
claw  probe  interference  effects. 


SECTION  VII 

SHOCK  ORIENTATION  EFFECTS 


k 


This  section  considers  briefly  some  gross  effects  of  the  blast 
intercept  angle,  <f>,  on  engine  face  total  pressures  for  the  test  model 
as  calculated  from  the  BID  code  for  Mach  0.85  flight  conditions. 

It  was  noted  in  Section  5 that  for  a blast  wave  of  constant  over- 
pressure, the  increase  in  blast  total  pressure  in  free  space  (before 
intercepting  the  inlet)  is  largest  for  zero  intercept  angle  and  decreases 
with  increasing  intercept  angle,  according  to  the  curves  of  Figure  5.2. 

As  might  be  expected,  the  same  type  of  variation  is  obtained  for  the 
total  pressure  at  the  engine  face,  as  is  shown  in  Figure  7.1.  This 
figure  presents  total-pressure  time  histories  near  the  center  of  the 
blastward  and  leeward  inlets  as  calculated  with  the  BID  code  for  a 
constant-strength  5-psi  overpressure  blast  wave  striking  the  inlets 
at  intercept  angles  of  90°,  105°  and  135°. 

It  may  be  noted  in  Figure  7.1  that,  in  addition  to  a decrease  of 
face  total  pressure  with  increasing  intercept  angle,  there  is  a noticeable 
reduction  in  the  complexity  of  the  pressure  signature  for  the  blastward 
inlet  as  the  intercept  angle  increases  from  90°  to  135°.  For  the  nearly 
side-on  intercept  angles  of  90°  and  105°  there  are  two  distinct  peaks, 
the  first  representing  the  initial  blast  wave  striking  the  face  and  the 
second  representing  the  reflected  shock  from  the  choked  section  of  the 
simulated  engine.  For  the  larger  135°  intercept  angle,  where  the  blast 
wave  overtakes  the  inlet  from  the  rear,  the  distinction  between  the 
incident  and  reflected  blast  wave  is  hardly  detectable. 

It  can  also  be  observed  in  Figure  7.1  that  the  rate  of  increase  of 
face  pressure  to  its  maximum  value  decreases  considerably  as  the  intercept 
angle  increases  from  90°  to  135°,  both  for  the  blastward  and  leeward  inlets. 

In  order  to  better  relate  the  face  pressures  to  the  incident  blast 
conditions,  the  total-pressure  curves  from  Figure  7.1  are  also  presented 
in  Figure  7.2  in  the  form  of  the  ratio  of  the  change  in  total  pressure 
(due  to  the  blast)  to  the  change  in  total  pressure  across  the  blast  front 
before  it  encounters  the  inlet  (from  Figure  5.2).  It  is  evident  from 

132 


__________ 


y' 


Comparison  of  engine  face  total  pressure  time  histories  for  three  blast  intercept  angl 


Comparison  of  engine  face  total  pressure  ratio  time  histories  for  three  blast 
intercept  angles. 


(b)  Leeward  inlet. 


Figure  7.2  that,  for  the  blastward  inlet,  the  average  total-pressure 
change  at  the  simulated  engine  face  has  about  the  same  ratio  to  the 
incident  total-pressure  change  for  all  three  intercept  angles,  the  ratio 
being  about  1.17  for  the  time  range  shown.  For  the  leeward  inlet,  this 
ratio  is  smaller  and  varies  somewhat  with  the  intercept  angle. 

The  significance  of  the  above  observations  with  respect  to  engine 
performance  has  not  yet  been  explored  in  detail.  However,  both  the 
decreasing  intensity  and  the  decreasing  rate  of  change  of  total  pressure 
observed  for  increasing  intercept  angle  in  Figure  7.1  may  imply  a 
decrease  in  inlet  distortion  with  increasing  intercept  angle,  at  least 
in  the  90°  to  135°  range. 
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SECTION  VIII 

INTERACTION  OF  BLAST  SHOCK  WITH  THE 
ENGINE  FAN 


The  principal  question  regarding  blast  wave  interaction  with  the 
engine  inlet,  of  course,  is  what  effect  this  interaction  would  have 
upon  the  operation  of  the  engine  downstream. 

The  effect  would  depend  upon  the  nature  of  the  transient  events 
that  take  place  when  the  blast  shock  and  blast-induced  flow  reach  the 
engine.  At  present  there  is  no  means  for  analyzing  this  sequence  of 
events  nor  are  there  test  data  available  for  its  assessment.  As  an 
admittedly  limited  first  step  in  understanding  the  problem,  the  inter- 
action at  the  fan  stage  of  a turbofan  engine  will  be  examined  here  on 
a quasi-steady  basis. 

8.1  QUASI-STEADY  BLAST  SHOCK  INTERACTION  WITH  A FAN 

The  quasi-steady  interaction  of  a blast  shock  leaving  the  inlet 
and  entering  a fan  will  be  analyzed  by  making  the  assumption  that  the 
interaction  reaches  an  equilibrium  rapidly  enough  that  a new  steady-state 
fan  point  is  reached  essentially  immediately  after  shock  arrival.  Of 
course,  the  transient  interactions  of  the  blast  shock  with  the  rotor 
and  stator  blades  are  expected  to  be  important,  but,  until  the  means 
are  available  for  analysis  of  these  interactions,  the  analysis  of  a 
quasi-steady  interaction  will  have  to  suffice  for  providing  an  assessment 
of  some  of  the  features  of  shock  interaction  with  a fan. 

The  time  sequence  of  events  associated  with  the  interaction  of  a 
shock  wave  with  a fan  is  shown  in  Figure  8.1.  Here  the  flow  going  into 
and  out  of  the  fan  is  treated  on  a one-dimensional  basis.  Region  2 
contains  the  preblast  flow  entering  the  fan  and  region  3 contains  the 
preblast  flow  leaving  it.  The  fan  is  located  at  the  origin  of  the 
coordinates . 


shock.  The  shock  wave  is  partially  reflected  upstream  producing  the 
flow  in  regions  2"  ahead  of  the  fan.  It  also  is  transmitted  through 
the  fan  producing  the  flow  in  regions  3'  and  3"  downstream  of  the  fan. 

The  air  in  region  3'  has  first  passed  through  the  fan  and  is  then 
shock  compressed  by  the  transmitted  shock.  The  air  in  region  3"  has 
been  twice  shock  compressed  before  passing  through  the  fan,  first  by 
the  incident  shock  (separating  zones  2 and  2')  and  then  by  the  reflected 
shock  (separating  zones  2'  and  2").  The  air  in  regions  3'  and  3"  is 
separated  by  a contact  surface;  the  pressures  and  particle  velocities 
are  equal  in  the  two  regions  but  the  densities  differ. 

It  should  be  noted  that  the  shock  interaction  with  a fan  might 
possibly  produce  expansion  waves  in  place  of  the  transmitted  or  reflected 
shock  waves.  For  the  particular  fan  interactions  analyzed  here,  which 
are  believed  to  be  typical  for  modern  axial-flow  turbo-fan  engines, 
only  shock  waves  would  be  possible,  so  only  shock  waves  are  assumed 
here  for  the  transmitted  and  reflected  waves.  The  addition  of  expansion 
waves  to  the  analysis  would  be  straightforward. 

Prior  to  the  arrival  of  the  incident  shock,  the  fan  compresses 
region-2  flow  upstream  into  region-3  flow  downstream.  After  arrival 
of  the  incident  shock  at  the  fan,  the  fan  compresses  region-2"  flow 
upstream  into  region-3"  flow  downstream.  The  fan  therefore  goes  from 
a pre-interaction  operating  point  between  regions  2 and  3 to  a post- 
interaction operating  point  between  regions  2"  and  3". 

The  flow  properties  in  the  initial  regions  2 and  3 are  determined 
by  a match  between  the  fan  characteristics  and  the  inlet  and  downstream 
engine  characteristics.  For  typical  current  aircraft  the  flows  entering 
and  leaving  the  fan  are  subsonic,  so  the  inlet,  fan  and  downstream 
engine  characteristics  are  interacting.  The  flow  properties  in  region 
2'  are  determined  by  knowing  the  strength  of  the  incident  shock.  The 
properties  in  regions  2",  3'  and  3"  are  determined  by  a match  between 
the  strengths  of  the  transmitted  and  reflected  shock  waves  and  the 
characteristics  of  the  fan.  Results  for  a typical  fan  will  be  shown. 


This  quasi-steady  interaction  of  the  shock  with  a fan  has  been 
programmed  into  a FORTRAN  IV  code  called  QBIF  (Quasi-steady  Blast 
Interaction  with  a Fan).  The  fan  characteristics  are  a code  input. 

The  important  questions  stemming  from  the  QBIF  calculation  are 
(a)  what  operating  point  the  fan  would  arrive  at,  (b)  what  effect  the 
reflected  shock  would  have  upon  upstream  flow  within  the  inlet  and 
(c)  what  effect  the  transmitted  shock  wave  might  have  on  the  flow 
downstream  at  the  compressor  inlet  or  through  the  bypass  duct  to  the 
afterburner  and  the  turbine  outlet. 

8.2  BLAST  SHOCK  AT  ENGINE  FACE 

In  view  of  the  large  number  of  tests  and  variables  in  the  0.1-scale 
B-l  inlet  tests  at  AEDC,  calculations  with  the  QBIF  code  have  been  con- 
centrated in  this  study  on  the  test  conditions  of  a selected  run  in 
the  test  series  at  AEDC,  Run  40  (Part  619).  The  pre-blast  flow  and  blast 
conditions  for  this  run  are  tabulated  in  Table  8.1. 

The  wind-tunnel  Mach  number  for  Run  40  (Part  619)  was  0.85  and  the 
mass  flow  at  the  engine  faces  was  351.9  and  351.0  lb/sec  full-scale 
for  the  outboard  and  inboard  inlets,  respectively.  The  blast  shock 
strength  scaled  to  standard  sea  level  conditions  was  5.8  psi.  The  jump 
in  average  total  pressure  at  the  two  engine  faces  when  the  blast  shock 
arrived  was  0.435  and  0.150  of  the  pre-blast  free-stream  total  pressure, 
respectively,  for  the  two  inlets.  The  corresponding  ratios  of  the 
static  pressure  behind  the  shock  to  the  value  ahead  is  1.280  and  1.101 
at  the  engine  face  for  the  respective  inlets. 

8.3  ENGINE  FAN  CHARACTERISTICS 

For  two-stage  fan  characteristics  representative  of  current 
technology,  the  data  of  Reference  8.1  have  been  employed.  A sketch  of 
this  fan  and  a tabulation  of  its  passage  dimensions  are  reproduced  from 
Reference  8.1  in  Figure  8.2.  The  fan  was  designed  for  a tip  speed  of 
1450  ft/sec,  an  overall  pressure  ratio  of  2.8  and  a corrected  mass  flow 
of  184.2  lbm/sec.  The  efficiency  at  this  design  operating  point  was 
found  to  be  85.7  percent  from  tests  (Reference  8.1).  The  overall 
performance  maps  for  the  fan  are  presented  in  Figure  8.3.  An  enlarged 
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TABLE  8.1 


INLET  DATA  OF  RUN  40  (PART  619)  PERTINENT  TO 
BLAST  SHOCK- FAN  INTERACTION 


TUNNEL  DATA 

PREBLAST 


M =0.850 
o 

Pt  = 11.801  psia 
o 

p = 7.358  psia 
ro 

T = 569°F 
o 


ENGINE-FACE  DATA 


OB  IB 

INLET  INLET 

PREBLAST 

M2  0.519  0.511 

p /p  0.979  0.980 

av  o 

WZR-FS (lb/sec)  351.9  351.0 

BLAST 

p -p  )/p  0.435  0.150 

av  av  o 
s 

pg/p  1.280  1.101 


Scaled  to  one  atmosphere  ambient  pressure. 


BLAST 


Ap  = 5.8  psi 
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(b)  Schematic  of  two-stage  fan  test  arrangement. 
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Figure  8.2.  Configuration  of  NASA  two-stage  axial-flow  fan. 
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Two  stage  fan  overall  performance 


Figure  8.3.  Performance  maps  for  NASA  two-stage  axial-flow  fan 


map  is  presented  in  Figure  8.4  for  the  fan  total-pressure  ratio.  Four 
operational  lines  have  been  added  to  this  map  for  fan  exit  Mach  numbers 
of  0.45,  0.51,  0.64  and  1.00. 

Six  preblast  operating  points  were  calculated  for  the  two  engine- 
face  Mach  numbers  of  Table  8.1  and  for  three  values  of  the  stall  margin 
— 8,  15  and  24  percent.  The  definition  of  stall  margin,  SM,  used  here  is 


where  p /p 

(pt3/pt2?stall 
mass  flow. 


is  the  fan  total-pressure 
is  the  value  at  stall  or 


ratio  at  the  operating  point  and 
surge  for  the  same  equivalent 


8.4  Blast  Effect  on  Fan  Operation 

Calculations  with  the  QBIF  code  were  carried  out  for  the  fan  of 
Reference  8.1  for  three  initial  conditions  of  the  two  inlets.  The  fan 
operating  points  that  would  result  are  presented  in  Figure  8.4  for  the 
three  initial  stall  margins  for  the  two  inlets.  The  shock  pressure 
ratio  at  the  engine  face,  p2/p2>  is  the  parameter  for  the  curves.  The 
results  show  that  increasing  this  shock  strength  reduces  both  the  equivalent 
mass  flow  and  the  total  pressure  ratio. 

The  incident  shock  pressure  ratio  at  the  engine  face  for  the  inboard 
inlet  in  Part  619  was  1.101.  The  operating  points  resulting  from  this 
interaction  with  the  fan  for  the  three  initial  stall  margins  are  joined 
in  Figure  8.4  by  the  curve  for  p2/p2  = 1*10.  The  points  for  double  and 
half  that  shock  overpressure,  i.e.  p2/p2  e9ual  1*20  and  1.05,  are  joined 
similarly.  These  resulting  operating  points  are  all  well  within  the 
area  of  the  performance  map  where  acceptable  operation  would  be  expected. 
Therefore  the  blast  should  produce  no  problem  for  the  inboard  fan 
according  to  the  QBIF  code  calculations. 
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The  pressure  ratio  of  the  shock  at  the  engine  face  in  the  outboard 
inlet  in  Run  40  was  1.28.  The  resulting  three  operating  points  for  the 
three  stall  margins  wind  up  much  closer  to  the  M^=1.0  operating  line. 

The  flow  in  stator  no.  2 of  the  fan  would  be  expected  to  choke  at  about 
midway  between  the  M^=0.64  line  and  the  M^=1.00  line.  Therefore  the 
blast  would  be  expected  to  choke  the  fan  in  the  outboard  inlet  for  this 
shock  strength  and  an  initial  stall  margin  of  24  percent. 

The  resulting  operating  points  are  also  shown  for  about  a 42  and 
100  percent  increase  in  shock  overpressure,  i.e.  p^/P2  = 1-40  and  1.56, 
and  one-half  the  overpressure,  ^'2^2  ~ The  higher  overpressures 

would  be  expected  to  choke  the  fan  for  the  initial  24  percent  stall  margin. 

Choking  of  a fan  for  the  short  durations  that  would  be  experienced 
during  penetration  of  a blast  wave  is  not  expected  to  be  harmful  to  its 
operation  on  this  quasi-steady  basis.  The  problem  for  the  fan  would 
arise  if  it  be  driven  to  the  surge-stall  boundary.  The  effect  of  this 
initial  shock  from  the  blast  wave  for  the  cases  run,  however,  is  to  in- 
crease the  downstream  Mach  number,  thereby  driving  the  fan  away  from  the 
surge-stall  boundary.  What  the  effect  would  be  for  other  initial  fan 
points  remains  to  be  explored. 

8.5  FAN  REFLECTED  SHOCK 

The  properties  of  the  shock  reflected  upstream  from  the  fan,  following 
the  Interaction  of  the  incident  shock  with  the  fan,  are  presented  in 
Figure  8.5  from  the  calculations  with  the  QBIF  code. 

In  Figure  8.5(a)  the  pressure  of  the  reflected  shock  in  excess  of 
the  initial  pressure  is  presented  as  a ratio  to  the  overpressure  of  the 
incident  shock,  V2  p2  . This  excess  pressure  ratio  is  plotted  as  a 

P2"P2 

function  of  the  pressure  ratio  of  the  incident  shock,  p’/p2*  A value  of 
the  ordinate  of  unity  corresponds  to  no  reflected  shock  and  a value  of 
two  indicates  an  overpressure  of  the  reflected  shock  equal  to  the  over- 
pressure of  the  incident  shock. 
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STALL  MARGIN 


INCIDENT  SHOCK  PRESSURE 


INCIDENT  SHOCK  PRESSURE  RATIO 


A common  method  for  simulating  an  engine  in  inlet  testing  is  to 
choke  the  flow  leaving  the  inlet  at  about  the  downstream  location  for 
the  engine.  This  is  done  by  using  a fixed  choked  orifice  or  nozzle  or 
adjustable  vanes.  The  latter  were  employed  with  the  0.1-scale  B-l  inlet 
and  were  located  near  the  fan  station. 

The  properties  of  shocks  reflected  upstream  from  a choked  orifice 
(or  nozzle  or  vanes)  and  fan  are  compared  in  Figure  8.5(a).  The  values 
of  the  incident  shock  pressure  ratio  at  the  engine  face  for  Run  40  are 
indicated  on  the  abscissa.  The  increment  in  excess  pressure  ratio  of 
the  reflected  shock  for  the  fan  is  markedly  greater  than  for  the  choked 
orifice.  Under  the  shock  conditions  of  the  outboard  inlet  in  Run  40  it 
is  33  to  45  percent  greater.  This  means  that  the  strength  of  the  shock 
reflected  upstream  from  the  fan  would  be  that  much  greater  with  a turbo- 
fan engine  than  was  experienced  in  the  tests  in  the  0.1-scale  B-l  inlet. 

The  partial  Mach  number,  M^' , of  the  flow  behind  (downstream)  the 
reflected  shock  and  entering  the  fan  or  a choked  orifice  is  shown  in 
Figure  8.5(b).  The  inflow  Mach  number  for  a choked  fan  would  be  constant. 
The  decrease  for  the  fan  is  due  primarily  to  the  increase  in  temperature 
and  speed  of  sound  by  the  double  shocking  of  the  air,  and  only  partially 
to  the  reduction,  generally,  in  particle  velocity. 

Analyses  of  aircraft  turbine  engines  are  generally  concerned  with 
the  total  pressure.  The  excess  total-pressure  ratio, 

„ll  „ 
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2 2 , is  presented  in  Figure  8.5  (c) . It  is  not  as  large  for 


Pt  " Pt 
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either  the  choked  orifice  or  fan  as  was  the  excess  static  pressure.  In 
other  words,  this  means  that  the  total  pressure  does  not  get  incremented 
as  much  by  the  reflected  shock  relative  to  the  increment  by  the  incident 
shock  as  does  the  static  pressure.  However,  the  total  pressure  increment 
due  to  the  reflected  shock  is  markedly  greater  for  the  fan  than  for  a 
choked  orifice.  For  the  outboard  inlet  and  the  shock  conditions  of 
Run  40  (Part  619) , the  increment  would  be  34  to  50  percent  greater  for 
the  fan  than  for  the  choked  orifice. 
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The  variation  of  the  properties  of  the  shock  transmitted  downstream 
by  the  fan  as  a function  of  the  incident  shock  strength  are  presented  in 
Figure  8.6.  The  pressure  ratio  for  the  transmitted  shock  wave  is  shown 
in  Figure  8.6(a).  The  ratio  is  always  less  than  for  the  incident  shock, 
but  it  is  relatively  large.  The  effect  upon  the  downstream  engine 
operation  might  be  significant. 

The  particle  Mach  number  M!^  downstream  of  the  fan  for  the  transmitted 
shock  is  plotted  in  Figure  8.6(b).  It  increases  with  the  strength  of 
the  incident  shock  in  all  cases.  It  is  reasonable  to  expect  that  the 
second  stator  of  the  fan  would  choke  at  a downstream  Mach  number  somewhere 
near  0.8,  so  the  downstream  Mach  number  would  not  increase  beyond  that 
point.  In  the  absence  of  data  on  the  choke  Mach  number,  the  curves  have 
been  continued  as  if  there  were  no  choking. 


This  condition  for  a choked  second  stator  would  be  nearly  reached 
for  the  fan  operating  at  a stall  margin  of  24  percent  and  incident  shock 
pressure  ratios  approaching  1.6.  For  lower  initial  stall  margins  the 
pressure  ratio  of  the  incident  shock  would  need  to  be  considerably 
greater  for  choking.  In  summary,  the  result  of  this  analysis  with  the 
QBIF  code  is  to  indicate  that  this  fan  operating  at  the  initial  conditions 
of  the  Run  40  test  would  not  experience  any  operational  problems  on  a 
quasi-steady  basis  for  (step)  shocks  at  the  engine  face  of  the  magnitude 
expected  of  Run  40  for  the  incidence  angle,  <j> » of  100  degrees  and  shock 
overpressure  of  5.8  psi  scaled  to  1 atm.  The  upstream  reflected  shock 
would  be  notedly  stronger,  however,  than  for  the  choked  control  vanes 
used  (as  the  simulated  engine)  in  the  tests. 


The  effect  of  the  reflected  shock  on  the  flow  within  the  inlet  will 
be  examined  in  Section  9.  The  strength  of  the  transmitted  wave  is  large. 
Its  effect  on  engine  operation  downstream  is  not  known. 
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SECTION  IX 
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REFLECTED  SHOCK-BOUNDARY  LAYER  INTERACTION  WITHIN  INLET 


Another  important  feature  of  the  inlet-blast  interaction  problem  is 
the  interaction  inside  the  inlet  between  the  blast  shock  waves  and  the 
inlet  boundary  layer.  The  sketch  in  Figure  9.1  depicts  the  features  of 
shock-boundary  layer  interaction  that  might  occur  on  the  walls  of  the 
inlets  of  the  B-l  type  due  to  the  reflected  shock  from  the  fan.  The 
flow  in  the  boundary  layers  would  be  forced  to  negotiate  the  adverse 
pressure  gradients  associated  with  these  shock  waves  moving  upstream 
from  the  engine  forces.  If  sufficiently  strong,  the  adverse  pressure 
gradients  could  cause  the  boundary  layers  to  separate  on  the  inlet 
walls,  depending  upon  the  boundary  layer  characteristics.  This  pos- 
sibility will  be  assessed  in  this  section. 

The  interaction  is  complicated  by  the  variation  of  the  "free- 
stream"  properties  within  the  inlet  both  in  the  flow  direction  and  time- 
wise.  The  total  pressure  and  total  temperature  in  particular  both  vary 
longitudinally  and  with  time  as  the  reflected  shock  moves  upstream. 

More  data  are  required  for  this  analysis  than  were  measured  in  the  tests 
with  the  0.1-scale  B-l  inlet,  so  calculations  made  of  the  flow  with  the 
BID  code  have  been  employed. 

9.1  SHOCK-BOUNDARY  LAYER  CALCULATION 


The  development  of  the  boundary  layers  was  computed  using  the 
results  of  BID  code  calculations  to  provide  the  distributions  of  the 
free-stream  properties  of  static  pressure,  total  pressure  and  total 
temperature  near  each  wall.  The  BID  code  models  the  fan  as  a choked 
orifice,  in  line  with  the  0.1-scale  inlet  model.  The  distributions  for 
each  wall  of  the  inlets  were  selected  corresponding  to  the  time  when  the 
static  pressure  gradients  appeared  to  be  relatively  steep  within  the 
blastward  inlet.  No  systematic  attempt  was  made  to  select  the  most 
adverse  distributions  for  the  boundary  layers.  These  properties  defined 
the  conditions  at  the  outer  edge  of  each  boundary  layer. 


I 
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The  growth  and  behavior  of  the  boundary  layers  along  the  walls  for 
these  free-stream  conditions  were  computed  using  the  BLAYER  code  of 
NASA,  presented  in  Reference  9.1.  Since  the  Reynolds  number  of  the 
boundary  layers  would  be  generally  large  relative  to  the  transition 
value,  the  boundary  layers  were  assumed  to  be  turbulent  from  the  leading 
edge  of  each  surface. 

The  BLAYER  code  assumes  a constant  total  pressure  and  total  tempera- 
ture outside  of  the  boundary  layer.  Because  they  both  generally  vary 
significantly  for  the  calculations  made,  the  calculations  were  carried 
out  piecewise  taking  the  total  pressure  and  total  temperature  as  constant 
over  each  segment.  This  procedure  was  found  to  be  significant  to  the 
values  computed,  but  it  is  not  believed  to  alter  the  general  conclusions 
reached  from  the  calculations. 

The  distributions  of  the  free-stream  properties  that  were  selected 
were  assumed  to  be  frozen.  That  is  they  were  taken  as  non-varying  in 
time,  whereas  they  do  change  as  the  reflected  shock  moves  upstream.  The 
significance  of  this  assumption  is  not  readily  assessable.  The  BLAYER 
code  applies  only  for  the  frozen  pressure  distribution,  so  the  question 
must  be  resolved  at  a later  time. 

The  BLAYER  calculations  were  performed  for  the  full-scale  B-l 
inlets,  as  scaled  up  from  the  0.1-scale  model.  The  blast  wave  in  the 
BID  calculation  is  assumed  to  arrive  from  the  side  of  the  aircraft,  that 
is  = 90°  (0°  is  a head-on  intercept),  with  a shock  overpressure  of  5.0 
psi.  The  flight  Mach  number  was  taken  to  be  0.85  and  the  mass  flow  rate 
in  the  inlets  to  be  350  lb/sec. 

9.2  SHOCK-BOUNDARY  LAYER  RESULTS 

The  pressure  distributions  for  the  boundary  layer  calculations  with 
the  BLAYER  code  are  presented  in  Figure  9.2  along  with  the  boundary- 
layer  incompressible  form  factor  that  resulted  from  these  calculations. 

The  boundary-layer  factor  H is  a parameter  essentially  representing 
the  velocity  distribution  across  the  boundary  layer  (normal  to  the  wall) . 
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When  the  velocity  distribution  is  more  uniform,  the  parameter  approaches 
unity.  A value  of  1.2  to  1.4  is  typical  for  a turbulent  boundary  layer 
on  a flat  plate  in  the  absence  of  a pressure  gradient.  In  an  adverse 
pressure  gradient  (increasing  pressure  in  the  direction  of  flow)  the 
velocities  decrease  toward  the  wall  so  increases.  When  the  velocity 
gradient  normal  to  the  wall  goes  to  zero,  the  flow  can  separate  from 
the  wall.  Experience  shows  that  separation  may  occur  for  values  as 
low  as  2.0  and  definitely  by  2.8. 

The  th  distributions  along  the  walls  for  the  blastward  inlet  are 
presented  in  Figure  9.2(a).  The  results  indicate  that  the  boundary 
layers  would  separate  on  both  the  cowl  and  the  ramp  surfaces.  The  cowl 
boundary  layer  would  separate  upstream  of  the  ramp,  so  it  is  quite 
probable  that  the  stream  would  actually  be  deflected  to  the  ramp  side 
and  that  the  ramp  boundary  layer  would  not  separate.  In  any  event, 
separation  would  take  place  on  one  side  of  the  blastward  inlet  which 
would  produce  significant  distortion  at  the  engine  face  when  the  separated 
flow  arrives  there. 

The  distribution  for  the  leeward  inlet  is  shown  in  Figure  9.2(b). 
The  distribution  for  the  ramp  is  not  presented  because  the  longitudinal 
variation  in  total  pressure  outside  of  the  boundary  layer  is  so  large  as 
to  make  conclusions  from  the  results  undependable.  The  distribution 
for  the  cowl  boundary  layer  reaches  a maximum  of  only  1.8,  which  is 
below  the  threshold  of  separation. 

In  the  BID  code  the  fan  is  modeled  as  a choked  orifice  as  mentioned 
(for  matching  to  the  0.1-scale  B-l  inlet  model  data),  so  pressure  gradients 
for  the  leeward  cowl  may  be  15  to  40  percent  greater  with  a fan.  This 
might  be  enough  to  separate  the  cowl  boundary  layer  but  the  question  is 
unresolved. 

There  is  an  effective  antificial  viscosity  in  the  BID  code  because 
the  cells  are  fixed  in  space  (also  called  "shock  capturing").  This  has 
the  effect  of  reducing  pressure  gradients,  resulting  in  lower  calculated 
values  of  H^. 
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The  question  then  arises  as  to  why  larger  values  of  the  distortion 
parameters  were  not  measured  during  the  test  periods  with  the  0.1-scale 
B-l  inlet.  Some  rough  estimates  of  the  time-distance  relationships  of 
events  within  an  inlet  are  useful  in  examining  this  question.  Several 
are  depicted  in  Figure  9.3  for  Mach  0.85  flight  conditions.  The  initial 
shock  would  reach  the  engine  face  in  about  0.9  milliseconds  and  the 
vanes  is  about  1.1  milliseconds.  The  reflected  shock  from  the  vanes 
moves  upstream  more  slowly,  as  shown.  Somewhere  within  the  inlet  the 
boundary  layer  would  separate.  To  estimate  when  the  separated  flow 
might  reach  the  engine  face,  some  other  velocities  must  be  considered. 

The  particle  velocity  within  the  duct  is  about  one-third  of  the 
initial  (incident)  shock  velocity.  Vorticity  shed  from  the  leading  edge 
of  the  cowl  or  splitter  lips  would  first  travel  at  about  one  half  of  the 
particle  speed  and  slowly  speed  up  as  it  moves  out  into  the  free  stream. 
Assuming  a step  speed-up  half  way  down  the  duct,  the  vortex  would  reach 
the  engine  face  at  4.2  milliseconds. 

The  blast-type  flow  period  produced  by  the  shock  tube  exit  flow  at  the 
inlet  is  terminated  by  the  arrival  of  the  shock-tube  cold  gas  or  the 
backward-facing  shock  or  compression  wave,  observed  as  a sudden  increase 
in  the  rate  of  pressure  change  at  the  inlet.  In  the  case  shown  in 
Figure  9.3  this  event  would  occur  at  the  inlet  entrance  at  about  3.3 
milliseconds,  and  could  arrive  at  the  engine  face  by  about  4.2  milliseconds. 

The  speed  of  the  front  of  the  separated  boundary  layer  is  not 
known,  but  it  is  reasonable  to  expect  it  to  be  significantly  less  than 
for  the  free-stream  particles.  For  estimation  purposes,  if  the  reflected 
shock  would  separate  the  boundary  layer  midway  down  the  inlet  and  the 
front  of  the  separated  flow  travels  at  one-half  of  the  particle  velocity, 
the  separation  zone  would  reach  the  engine  face  at  somewhat  more  than  a 
millisecond  after  the  termination  of  the  test  period. 

Figure  5.8  shows  that  the  distortion  parameter  IDC  increased  con- 
siderably for  the  blastward  inlet  for  Ap  = 5.0  and  5.2  psi  beginning  at 
about  one  millisecond  after  termination  of  the  3.3-millisecond  nominal 
test  period,  indicated  in  the  figure  by  the  arrows.  IDR  increased  at 
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Inlet  time-distance  relationships  for  a flight  Mach  number  of  0.85 
and  a mass  flow  rate  of  350  lb/sec. 


about  the  same  time  but  not  in  as  great  a proportion,  in  comparison  with 
the  values  during  the  test  period.  IDL  increased  apparently  well  beyond 
the  full-scale  value  of  1.50. 

For  the  leeward  inlet,  Figure  5.9,  the  distortion  indicated  is 
not  as  large.  IDC  and  IDL  appear  to  have  increased  beginning  at  about 
the  end  of  the  nominal  test  period,  particularly  for  the  5.2  psi  shock. 
IDR,  on  the  other  hand,  decreased.  IDL  experiences  only  brief  bursts 
above  a value  of  about  1.1. 

Because  these  large  increases  in  distortion  occurred  after  the 
expected  arrival  of  the  backward  facing  shock  or  compression  wave,  a 
question  does  arise  as  to  what  degree  the  high  level  of  distortion  may 
have  been  due  to  the  character  of  the  blast  wave  during  the  later  period. 
This  question  should  be  examined  to  the  extent  possible  for  the  present 
data. 


SECTION  X 
DISCUSSION 


This  section  presents  a summary  and  discussion  of  some  of  the 
important  aspects  of  the  development  and  testing  program,  most  of  which 
are  described  in  detail  in  previous  sections,  and  indicates  some  recommen- 
dations for  future  developments.  Topics  discussed  are  the  blast  generator 
development,  the  16T  wind  tunnel  blast  tests,  the  BID  code  for  predicting 
inlet  blast  flows  and  blast-fan  interactions. 

10.1  BLAST  GENERATORS 

The  blast  generation  development  program  began  with  theoretical 
REFLECT2/S2D  code  calculations  for  the  firing  of  a shock  tube  into  still 
air,  followed  by  firings  of  two-in. -dia  shock  tubes  into  the  AEDC  IT 
(one-foot-square)  wind  tunnel  at  Mach  numbers  from  0.6  to  0.9  (Sec.  2). 
Results  of  these  tests,  for  the  case  of  the  shock  tube  firing  perpendicular 
to  the  tunnel  axis,  gave  useful  blast  type  flows  in  this  Mach  number 
range  for  the  overpressure  range  2-5  psi  (referred  to  sea-level  ambient 
pressure)  for  distances  from  the  tube  exit  of  about  4.5  diameters,  for 
polar  angles  from  the  tube  axis  from  -15°  to  +30°  downstream  of  the 
tube.  Durations  of  nearly  constant  blast  pressure  (to  30%  decay)  were 
obtained  for  these  conditions  in  the  range  of  about  0.5  to  2 milliseconds 
per  foot  of  tube  diameter,  with  best  results  being  obtained  at  a polar 
angle  of  15°.  At  downstream  polar  angles  from  the  tubes  of  45°  or  more, 
the  blast  wave  durations  were  generally  much  shorter  and  were  not  of 
interest.  Also,  it  was  found  from  the  IT  tests  that  if  the  shock  tube 
driver  pressure  was  too  low  compared  to  the  ambient  pressure  a non- 
constant pressure  type  of  blast  wave  profile  was  produced,  characterized 
by  an  extremely  rapid  decay. 

Results  of  the  above  calculations  and  small  shock  tubes  were  found 
to  be  essential  in  selecting  the  shock-tube  dimensions,  locations  and 
driver  pressures  for  the  subsequent  development  of  large  shock  tubes 
for  the  AEDC  16T  (16  foot  square)  wind  tunnel  tests. 

The  three  large  shock  tube  installation  designed  and  constructed 
for  the  16T  wind  tunnel  tests  had  the  following  characteristics. 

Distances  from  the  shock  tube  exits  to  the  inlet  were  3.8  to  4.3  tube 
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diameters,  at  polar  angles  from  11.2°  to  28.3°.  Shock  tube  driver 
pressures  of  55  to  186  psia  were  used,  usually  at  a tunnel  ambient 
pressure  of  one-half  atmosphere. 

Definition  of  the  properties  of  the  blast  waves  produced  at  the 
inlet  location  in  the  16T  tests  was  made  using  both  measured  data  from 
the  tests  as  well  as  the  test  results  from  the  IT  tests.  The  strength 
of  the  blast  shock  was  well  defined  at  probe  locations  by  the  16T 
measurements,  but  interference  from  the  inlet  model  degraded  the  definition 
of  the  flow  following  the  shock  somewhat.  Firings  in  the  16T  tunnel 
without  a model  would  be  useful  for  defining  these  blast  properties, 
but  they  were  precluded  by  limitations  on  tunnel  time;  as  an  alternative, 
an  analysis  of  model  interference  on  the  blast  measurements  is  recommended. 
The  incidence  angle  of  the  blast  shock  was  defined  from  comparative 
arrival  times  at  several  pressure  transducers;  shadowgraphs  would  be  a 
useful  alternative  for  this  measurement  in  future  tests  for  greater 
accuracy. 

Shock-tube  dimensions  of  a 22.6  in-ID  and  a 17-ft  length  were 
selected  to  provide  the  best  combination  of  blast  wave  characteristics 
within  space  constraints  of  the  wind  tunnel  for  the  three  shock  tubes. 
Commercial  pre-scored,  pre-stressed  diaphrapms  (rupture  discs)  in  a 
double-diaphragm  arrangement  provided  fast,  clean-opening  diaphragm 
bursts  free  of  any  petal  debris.  By  using  three  shock  tubes  the  number 
of  firings  during  the  test  period  was  essentially  doubled  over  using  a 
single  tube. 

10.2  THE  16T  WIND  TUNNEL  TESTS  AND  ANALYSIS 

The  0.1-scale  B-l  inlet  model  used  in  the  16T  tests  had  been  mated 
previously  to  the  tunnel  sting  support,  so  only  some  strengthening  of 
the  model  and  the  16T  tunnel  walls  for  loading  by  the  shock-tube  driver 
gas  was  required  for  testing. 

Blast  tests  were  performed  at  Mach  0.55,  0.70,  0.85  and  0.90  with 
inlet  mass  flow  rates  ranging  from  235  to  350  lb/sec  (full  scale), 
shock  overpressures  from  about  2 to  6 psi  (referred  to  sea  level)  and 
model  yaw  angles  of  0 and  5 degrees  (nose  away  from  shock  tube) . 
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Typically,  blast  type  flows  lasting  about  3.3  milliseconds  were 
achieved  at  the  inlet,  with  corresponding  nearly  constant  pressure 
durations  (to  30%  decay)  of  between  1 and  4 msec.  This  met  test 
requirements  amply.  However,  for  future  tests  vortex  and  boundary-layer 
separation  phenomena  effects  that  are  now  believed  to  be  important  at 
later  times  make  it  desirable  to  double  the  test  durations. 

The  test  durations  might  be  increased  about  50  percent  by  moving 
the  model  about  four  feet  off  the  centerline  of  tunnel,  increasing  the 
distance  from  the  shock  tubes  and  thereby  the  blast  duration,  but  the 
shock-tube  dimensions  would  need  to  be  reexamined.  Any  additional  increase 
in  duration  would  require  using  a smaller-scale  model  (gaining  thereby 
in  scaled  duration) . 

The  blast  interaction  with  the  inlet  as  observed  in  the  tests  can 
be  divided  into  three  characteristic  phases  in  terms  of  events  at  the 
engine  face.  During  the  first  phase  the  blast  shock  travels  down  both 
inelts  and  reflects  from  the  choked  vanes  (or  fan  for  the  full-scale 
engine),  then  returning  to  the  engine  face.  In  the  second  phase  any 
vortex  shed  from  the  leading  edge  of  the  splitter,  for  example,  would 
reach  the  engine  face,  and  shock  wave-boundary  layer  interaction  within 
the  inlet  would  take  place  and  also  reach  the  engine  face.  After  these 
two  stages,  events  are  believed  to  happen  relatively  slowly,  that  is 
on  a quasi-steady  basis,  so  steady-state  test  results  and  analyses  would 

be  applicable  this  is  called  the  third  stage.  Attention  here  is 

focused  primarily  on  the  first  two  stages. 

In  the  blastward  inlet  during  the  first  period  the  multiple  reflec- 
tions of  the  shock  result  in  a "staircase"  of  shocks  arriving  at  the 
engine  face  over  about  a 0.3-msec  period  followed  by  some  decay  preceding 
the  arrival  of  the  reflected  shock  in  about  another  millisecond.  In 
the  leeward  inlet  the  shock  diffracts  around  the  splitter  and  enters 
the  inlet  more  normal  to  the  model  centerline.  It  reaches  the  engine 
face  as  essentially  a single  shock,  much  weaker  than  for  the  blastward 
inlet . 
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When  the  shock  in  the  blastward  inlet  reaches  the  engine  face, 
distortion  IDL  values  near  unity  are  observed  which  might  be  of  some 
concern  during  maneuvering  or  other  stressful  flight  conditions.  The 
distortion  patterns  examined  have  had  high  total  pressures  on  the  cowl 
side  of  the  face.  The  distortions  in  the  leeward  inlet  during  this 
period  were  generally  lower. 

Calculations  were  made  of  the  blast  shock  interaction  with  a fan 
typical  of  current  technology  for  turbofan  engines.  Shock  strengths 
measured  at  the  engine  face  in  the  16T  tests  were  used.  The  interaction 
was  assumed  for  simplification  purposes  to  be  quasi-steady  with  a step 
shock  (no  ramp  or  decay) . For  the  blastward  inlet  the  results  indicate 
the  fan  would  choke  at  the  second  stator  for  the  lateral  (4>=“90°)  blast 
intercept.  This  presumably  presents  no  operational  problem,  but  the 
results  are  not  believed  to  be  general.  The  interaction  should  be  studied 
for  the  range  of  conditions  expected.  Of  clear  concern  to  engine  operation 
are  effects  produced  downstream  by  the  transmitted  shock  wave  upon  the 
afterburner  and  turbine  rear  stages  and  upstream  by  the  reflected  shock 
upon  flow  within  the  inlet. 

Concerning  the  latter,  boundary-layer  calculations  indicate  that  the 
reflected  shock  would  definitely  separate  the  upstream  boundary  layers 
within  the  blastward  inlet  for  a 5-psi  blast  shock.  This  is  expected  to 
result  in  serious  distortion  at  the  engine  face. 

In  the  tests  IDL  values  over  1.25  were  observed  at  the  engine  face 
at  about  4 milliseconds  after  shock  arrival  there.  This  is  about  one 
millisecond  after  the  end  of  the  validated  duration  of  blast  type  flow, 
therefore  it  is  possible  the  distortion  at  that  late  time  was  due  to 
variations  in  the  blast  wave.  In  view  of  the  significance  of  this  effect 
to  engine  operation,  it  is  recommended  that  the  blast  data  be  analyzed 
further  to  verify  that  the  late-time  distortions  are  relevant  to  the 
blast-inlet  problem. 

Interaction  of  the  reflected  shock  with  the  boundary  layers  in  the 
leeward  inlet  could  only  be  analyzed  to  a limited  extent,  because  of 
restrictions  in  the  methods  of  analysis  available.  Separation  on  the 
cowl  side  might  occur  if  calculations  were  refined,  but  not  as  carried 
out.  Estimation  of  separation  on  the  ramp  side  requires  an  assessment 
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of  separation  effects  in  BID  calculations  for  the  leading  edge  of  the 
ramp  and  the  improvement  in  the  boundary-layer  method,  needed  also  for 
the  blastward  inlet. 

Vortex  formation  at  the  splitter  is  another  possible  source  of 
distortion  for  the  leeward  inlet.  Distortion  values  in  the  leeward  inlet 
did  reach  large  values  in  the  tests  for  times  after  the  3.3-msec  blast 
flow  period  but  whether  the  shedding  of  vortices  is  also  a blast  inlet 
problem  depends  on  resolution  of  the  late-time  blast-wave  data. 

During  the  about  3.3-msec  period  of  blast  type  flow  the  maximum  IDL 
values  were  affected  by  the  inlet  mass  flow  rate.  But  there  were  no 
definite  significant  effects  of  either  Mach  number  or  inlet  yaw  angle 
on  IDL. 

The  whole  question  of  blast  distortion  relative  to  its  effect  upon 
engine  operation  is  an  open  question  at  this  point.  The  particular 
distortion  indices  used  in  work  with  the  B-l  inlet  (IDC,  IDR,  IDL  and 
IDT)  have  been  employed  here,  basically  as  parameters  that  have  met 
with  success  for  distortions  of  the  type  obtained  in  typical  flight 
situations  (maneuvering  flight,  throttle  transients,  etc.).  There  are 
many  other  distortion  parameters  in  use.  Also,  there  are  factors 
involved  in  a blast  interaction  that  are  not  reflected  in  these  parameters, 
such  as  the  rapidity  of  the  changes  due  to  the  shocks,  and  features 
of  the  spatial  distortions  that  are  significantly  different,  etc.  The 
effects  of  these  differences  on  engine  operation  must  be  answered  before 
the  question  can  be  resolved. 


The  BID  code  provides  good  predictions  of  the  observed  pressure 
histories  in  the  inlet  duct  and  at  the  engine  face,  particularly  in 
reproducing  the  principal  events,  such  as  shock  reflections.  The 
pressures  predicted  at  the  engine  face  tend  to  be  slightly  low  but  the 
trends  are  generally  good,  so  that  a simple  factor  (accounting  for 
numerical  viscosity  effects)  would  be  expected  to  bring  them  into  good 
agreement.  The  BID  code  at  present  does  not  include  blast  decay;  this 
addition  to  the  code  would  provide  an  improvement  at  late  times.  A 
second  improvement  that  is  recommended  is  to  extend  the  code  to  three 
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dimensions;  this  would  allow  for  modelling  the  hub  and  would  account  for 
differences  in  the  vertical  direction.  Work  is  recommended  to  further 
correlate  distortion  predictions  of  the  BID  code  with  the  test  data. 


BID  code  calculations  have  been  made  for  intercept  angles  from  90 
to  135  degrees  from  head  on.  The  90-deg  intercept  results  in  much  larger 
total-pressure  increment  at  the  engine  face  with  greater  pressure-time 
variations  than  for  the  more  rearward  intercepts.  Further  studies 
extending  the  blast  intercept  angles  are  recommended. 

Regarding  the  effects  of  the  observed  inlet  pressures  on  an  engine 
fan  stage,  the  calculation  of  the  quasi-steady  interaction  of  the  blast 
shock  with  a fan  and  with  a choked  flow  for  a 90-deg  intercept  indicates 
the  reflected  shock  would  be  considerably  stronger  for  a fan  then  for  a 
choked-flow  condition,  as  used  in  these  tests,  by  as  much  as  43  percent 
for  a typical  case.  Control  of  flow  by  means  of  a choked  nozzle  or  vanes 
is  a standard  technique  used  in  inlet  testing.  Consideration  should  be 


given  to  this  problem  of  weaker  reflected  shocks  in  planning  and  analyzing 


f 


SECTION  XI 
CONCLUSIONS 


From  the  analyses  and  wind-tunnel  tests  carried  out  to  simulate 
nuclear  blast  interaction  on  a B-l  type  engine  inlet,  the  following 
conclusions  are  reached. 


1.  Blast  simulations  capability  has  been  developed  and  demonstrated 
in  the  AEDC  16T  wind  tunnel  for  blast  strengths  of  2 to  6 psi, 
scaled  to  sea-level  conditions,  and  nearly  constant  blast 
pressures  (with  less  than  + 30-percent  decay)  lasting  from  1 to  4 
milliseconds . 

2.  A data  bank  of  blast  interaction  data  has  been  developed  with 
the  0.1-scale  B-l  inlet  for  lateral  intercepts,  varying 
wind-tunnel  Mach  number  from  0.55  to  0.90,  inlet  mass  flow 
rate  from  235  to  350  lb/sec.  (full  scale),  shock  overpressures 
from  2 to  6 psi  (sea  level)  and  model  yaw  angles  from  0 to 

5 degrees. 

3.  The  BID  code  predicts  satisfactorily  the  principal  features 
of  the  pressure  time  history  within  the  inlet  duct  and  at  the 
engine  face.  Additional  work  is  recommended  for  prediction  - of 
distortion  at  the  engine  face.  Recommended  additions  to  the 
code  include  accounting  for  blast  decay  and  extending  the  code 
to  three  dimensions. 

4.  Four  potentially  adverse  effects  to  engine  operation  from  blast 
interaction  were  identified:  blast-induced  distortion,  fan 
choking,  afterburner  blow  out  and  shock-boundary  layer  induced 
distortion. 

5.  Blast-induced  distortions  during  the  time  of  definite  blast-type 
flow  in  the  inlet  were  usually  smaller  than  but  did  sometimes 
reach  normal  inlet  allocation  levels  for  the  B-l  airplane. 
Distortions  at  later  times  greatly  exceeded  the  allocation. 
Boundary-layer  separation  in  the  inlet  by  simulated  fan  (vane) 
reflection  of  the  blast  shock  is  suspect  for  the  latter 
distortion  levels. 
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6.  Distortion  studies  made  here  are  severely  limited  by  the 

unavailability  of  a parameter  for  judging  the  effect  of  transient 
distortions  on  gas-turbine  engine  operation.  Blast  tests  with 
an  inlet-engine  combination  are  strongly  recommended  for  meeting 
this  need. 
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APPENDIX  A 

THEORETICAL  STUDIES  OF  THE  FLOW  FIELD  PRODUCED  BY 
FIRING  A SHOCK  TUBE  INTO  A STATIONARY  FLUID 

A. 1 INTRODUCTION 

In  order  to  obtain  background  and  planning  information  for  the  wind 
tunnel  tests  of  this  report,  KA  performed  a series  of  theoretical  studies 
for  the  axi-symmetrical  problem  where  a prescribed  flow  exits  from  a 
shock  tube  into  initially  stationary  air  (y  = 1.4)  at  an  ambient  pressure 
of  one  atmosphere  (14.7  psia) . The  external  transient  pressure  disturbance 
or  blast  field  produced  in  the  fluid  was  calculated  by  two  KA  computer 
codes  for  a wide  range  of  shock  tube  firing  conditions,  as  described  below. 

A.  2 S2D  CODE  STUDIES 


Initial  calculations  were  made  with  the  KA  S2D  Eulerian  computer 
code  (References  A.l  and  A. 2),  with  several  different  representations  of 
the  geometry  of  the  exit  end  of  the  tube.  In  one  representation,  shown 
in  Figure  A. la,  the  shock  tube  was  assumed  to  be  fired  into  a half-space 
bounded  by  a rigid  wall.  In  a second  representation,  the  shock  tube 
was  assumed  to  protrude  into  the  flow  from  a boundary  by  various  distances, 
X,  as  shown  in  Figure  A. lb,  with  a non-reflecting  boundary  condition 
being  assumed  to  apply  at  this  "boundary"  distance  behind  the  end  of  the 
shock  tube.  More  specifically,  the  flow  field  for  x < -X  in  Figure  A. lb 
was  assumed  to  have  properties  independent  of  x,  with  X being  varied  from 


about  0.3  to  1.5  tube  diameters. 

These  particular  S2D  solutions  were  used  to  determine  the  sensitivity 
of  the  external  flow  field  to  the  boundary  conditions  as  described  above. 
In  general,  it  was  found  that,  for  the  range  of  exit  pressures  considered, 
from  37  to  103  psia,  the  pressures  at  all  distances  two  or  more  tube 


diameters  downstream  of  the  tube  exit  were  essentially  independent  of 
the  assumed  boundary  conditions. 


Since  our  interest  was  primarily  in  large  downstream  distances,  the 
above  observations  indicated  that  we  could  more  profitably  study  this 
shock-tube  fluid  flow  problem  by  making  use  of  the  REFLECT2  computer  code 
described  below,  which  permits  better  resolution  of  pressure  transients. 


SHOCK  FRONT 


(a)  Wall  boundary. 


(b)  No  wall  boundary. 

r Mem  of  a shock  tube  firing  into  a stationary  fluid. 
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but  is  presently  less  capable  than  the  S2D  code  for  handling  complex 
boundary  conditions. 

A.  3 REFLECT2  CODE  STUDIES 

Calculations  were  made  with  the  KA  REFLECT2  Lagrangian  computer 
code  (References  A. 3 and  A. 4)  for  the  shock  tube  exit  flow  conditions  of 
Figure  A. la  for  a range  of  exit  pressures  (p£)  from  24  to  103  psia. 

Geometrical  properties  of  the  shock  tube  were  taken  as  : inside 
diameter  of  2 inches,  driver  length  (LR)  of  4.5  inches  and  driven  length 
(LN)  of  9 inches. 

The  shock  tube  exit  flow  conditions  of  pressure,  density  and  velocity 
were  taken  to  be  the  same  as  for  a perfect  one-dimensional  shock  tube 
flow  with  a specific  heat  ratio  of  y = 1*4.  In  some  solutions  the  exit 
flow  conditions  were  assumed  to  be  constant  corresponding  to  the  initial 
blast  field  behind  the  shock  front.  In  other  solutions,  account  was 
taken  of  the  change  in  density  produced  when  the  contact  surface  originally 
at  the  tube  diaphram  position  reached  the  exit  end  of  the  shock  tube. 
Consideration  of  this  contact  surface  change  appeared  to  have  little 
effect  on  the  initial  blast  phenomena  discussed  below  and  is  not 
discussed  herein. 

Figure  A. 2 presents  sample  time  histories  of  pressure  time  histories 
at  several  positions  along  the  axis  of  symmetry  (r  = radial  distance  from 
tube  exit,  d = diameter)  for  several  firing  pressures  (p^=  P^/po,P2i 
= ?2^0'  P4=  <*r^ver  Pressure>  = sh°ck  (exit)  pressure,  pq  = ambient 
pressure) . 

The  following  important  features  may  be  seen  in  Figure  A. 2.  For 
the  lowest  shown  driver  pressure  ratio,  p^  = 2.8,  the  blast  overpressure 
decays  very  rapidly  after  shock  arrival  reaching  zero  in  less  than  0.2  msec 
for  the  two  inch-diameter  shock  tube.  This  waveform  has  a too  rapid  decay 
to  be  of  interest  for  the  present  program.  For  the  intermediate  pressure 
ratio,  p^  = 10.4,  a much  longer  positive  overpressure  duation  is  obtained, 
about  0.4  msec  and  the  overpressure  is  much  more  nearly  constant  during 
the  positive  phase  than  for  p^  = 2.8.  Finally,  for  the  largest  pressure 
ratio,  p^  = 15.3,  the  positive  phase  duration  is  again  about  0.4  msec, 
and  the  overpressure  remains  nearly  constant  for  over  0.3  msec. 


Figure  A. 3 presents  the  calculated  initial  peak  or  transient  blast 
overpressure,  Ap,  in  the  fluid  along  the  axis  of  the  shock  tube,  as  a 
function  of  driver  pressure  ratio  p^,  and  exit  pressure  ratio,  p and 
radial  distance  in  tube  diameters  from  the  tube  exit,  r/d. 

The  initial  blast  pressure  is  observed  to  increase  nonlinearly  with 
increasing  pressure  ratio  and  to  vary  roughly  inversely  with  distance 
from  the  tube  exit  in  the  range  of  interest  (0.1  < Ap/p^  < 0.4). 

These  theoretical  results  can  be  represented  fairly  well  by  the 
empirical  equation: 

Ap/po  = 1.5  (6.39(r/d)p21_1'75  + 0.4)"1*5  (A.l) 

for  the  range  0.1  <Ap/pQ<1.0  where  l<r/d<10  and  2<p2^<7. 

A. 4 CONCLUSIONS 

The  above  results  of  these  studies  indicated  that  for  reasonable 
shock  tube  driving  pressure  (under  600  psia)  external  shock  overpressures 
in  the  desired  range  of  2 - 5 psi  could  be  obtained  at  distances  from  the 
tube  exit  in  the  range  of  3 - 10  diameters.  Also  the  calculated  durations 
of  the  shock  overpressure  appeared  to  be  long  enough  to  be  useful  for 
simulating  a nuclear  blast  on  a small  inlet  model.  For  example,  at  a 
driver  pressure  of  225  psia,  at  4.5  diameters  from  the  shock  tube  along 
the  tube  axis,  the  shock  overpressure  along  the  axis  of  the  tube  was 
calculated  to  be  about  3.0  psi  and  the  duration  of  nearly  constant 
overpressure  to  be  about  2.5  milliseconds  per  foot  of  tube  diameter. 


APPENDIX  B 

STEADY-STATE  PRE-BLAST  TEST  CONDITIONS  IN  THE  AEDC  16T  TUNNEL 


This  appendix  presents  AEDC  tabulations  of  the  steady-state  wind 
tunnel  and  inlet  operating  conditions  and  inlet  pressures  and  distortions 
existing  immediately  before  the  firing  of  each  shock  tube.  The  data  are 
presented  in  the  chronological  order  of  the  firings,  by  Part  number, 
and  are  preceded  by  a table  of  nomenclature. 

In  these  tables  the  following  parameters  are  in  error:  IDCO,  IDCI, 
IDLO  and  IDLI  for  data  set  2 and  PDO,  PDI,  PDPO  and  PDPI  for  both  data 
sets  1 and  2.  Corrected  values  of  these  parameters  are  available,  but  are 
not  included  here  since  they  are  not  essential  for  the  present  study. 
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Parameter 

AB14>-I 

AB2<J>-I 

ALPHA 

BETA (SIC) 

B5-1 

C 

C3 

CP1SS 

CP2SS 

CP3SS 

Date 

Day 

ERCODE 

E10 

G34-2 

H1I 

HL<J> 

H16-1 

HR 

I 


Steady-State  Printout  Nomenclature 
DNA  B-l  Blast  Effects  Test 
P41T-D2A  (TF-419) 

Description 

Outboard  and  inboard  total  bleed  zone  I areas, 

in. 2 

Outboard  and  inboard  total  bleed  zone  II 
areas,  in. 2 

Model  angle  of  attack  defined  positive 
nose  up,  deg 

The  angle  of  yaw  defined  positive  pilot's 
right,  deg 

Fuselage  configuration 

Structural  mode  control  fin's  angle  relative 
to  the  model  axis 

Structural  mode  control  fin  configuration 

Claw  probe  one  static  pressure  ratioed  to 
free-stream  total  pressure 

Claw  probe  two  static  pressure  ratioed  to 
free-stream  total  pressure 

Claw  probe  three  static  pressure  ratioed  to 
free-stream  total  pressure 

The  date  the  data  were  taken 

The  day  the  data  were  taken,  deg 

Tunnel  error  code 

Boundary  layer  control  system  configuration 
Wing  gutter  configuration 
Inboard  lip  height,  in. 

Outboard  lip  height,  in. 

Wind  hood  configuration 

The  hour  the  data  were  taken,  hr 
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Steady-State  Printout  Nomenclature  (Continued) 


Page  2 


Parameter 
IDA 
I DC 

IDL 

IDR 

IDT 

140 

IDC12 

IDC45 

J27 

L7 

M 

MA 

MB 

MIN 

MFR2 

MFR100 

MODE 

M2 

M2S 

M6 

NCN 

NSP 

N30-1 


i 

L 


Description 

Average  engine-face  distortion  index 


Engine-face  total  circumferential  distortion 
index 


Fan  stall  margin  ratio 

Engine-face  total  radial  distortion  index 
Engine-face  overall  total  pressure  distortion 
Inlet  configuration 


Engine-face  total  circumferential  distortion 
index  at  the  hub 


Engine-face  total  circumferential  distortion 
index  at  the  tip 


Porous  bleed  hole  configuration 

Cowl  configuration 

Free-stream  Mach  number 

First  ramp  Mach  number 

Mach  number  constants  box 

The  minute  the  data  were  taken,  min. 

Engine-face  mass  flow  ratio 


Engine-face  mass  flow  ratio  corrected  to 
100%  recovery 


Mode  code 


Engine- face  Mach  number  based  on  engine-face 
static  to  total  pressure  ratio 


Engine-face  Mach  number  based  on  weight  flow 

Flow  meter  configuration 

Nozzle  contour  number 

Normal  shock  parameter 

Nacelle  configuration 
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Steady-State  Printout  Nomenclature  (continued)  Page  3 


Parameter 

Description 

P 

Tunnel  static  pressure,  psfa 

Part 

Part  number 

PCA-1 

Primary  plenum  pressure,  psfa 

PCB-2 

Back-up  plenum  pressure,  psfa 

PD 

Back-up  prediction  of  peak  instantaneous 
fan  stall  margin  ratio 

PDP 

Primary  prediction  of  peak  instantaneous 
fan  stall  margin  ratio 

PLX 

Model  measured  pressure  ratioed  to  free-stream 
total  pressure  where  L refers  to  outboard 
( L= 1 ) and  inboard  (L=2)  and  X refers  to  the 
pressure  designation  (X=020,  030,  035,  etc.) 

PMIN 

The  minimum  engine-face  total  pressure 
ratioed  to  free-stream  total  pressure 

PMAX 

The  maximum  engine-face  total  pressure 
ratioed  to  free-stream  total  pressure 

PTI 

Tunnel  compressor  inlet  total  pressure,  psfa 

PTA-1 

Primary  tunnel  total  pressure,  psfa 

PTB-2 

Back-up  tunnel  total  pressure,  psfa 

PSTCEl 

Exit  control  pressure  for  shock  tube  number 
one,  psia 

PSTCE2 

Exit  control  pressure  for  shock  tube  number 
two,  psia 

PSTCE3 

Exit  control  pressure  for  shock  tube  number 
three,  psia 

PSTCDl 

Driver  control  pressure  for  shock  tube  number 
one,  psia 

PSTCD2 

Driver  control  pressure  for  shock  tube  number 
two,  psia 

PSTCD3 

Driver  control  pressure  for  shock  tube  number 
three,  psia 

PSTCIl 

Interstage  control  pressure  for  shock  tube 
number  one,  psia 
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Parameter 

PSTCI2 

PSTCI3 

PSTRD1 

PSTRD2 

PSTRD3 

PSTRI1 

PSTRI2 

PSTRI3 

P1700 

P1701 

Q 

R2 

RB<)>-I 
RC  «J>  — I 
RE  X 1 0 ~ 6 

SCHED 

Sec 

SHOCK  TUBE 

S19 

TAVG 

Test 

TI2 


Description 

Interstage  control  pressure  for  shock  tube 
number  two,  psia 

Interstage  control  pressure  for  shock  tube 
number  three,  psia 

Driver  pressure  for  shock  tube  one,  psia 

Driver  pressure  for  shock  tube  two,  psia 

Driver  pressure  for  shock  tube  three,  psia 

Interstage  pressure  for  shock  tube  one,  psia 

Interstage  pressure  for  shock  tube  two,  psia 

Interstage  pressure  for  shock  tube  three,  psia 

MA  total  pressure  ratioed  to  free-stream  total 
pressure 

MA  static  pressure  ratioed  to  free-stream  total 
pressure 

Free-stream  dynamic  pressure 

Engine-face  total  pressure  recovery 

Outboard  and  inboard  second  ramp  angle,  deg 

Outboard  and  inboard  third  ramp  angle,  deg 

Free-stream  unit  Reynolds  number  times  ten 
to  the  minus  six,  ft~l 

Tunnel  parameter 

The  second  the  data  were  taken,  sec 

SELECT  Shock  tube  select  code  with  the  last  two 
digits  referring  to  the  volumetric  size, 
cu. ft. 

Sideplate  configuration 

Engine-face  tip  average  static  pressure 
ratioed  to  free-stream  total  pressure 

Test  number 

Average  engine-face  turbulence  factor 
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Steady-State  Printout  Nomenclature  (concluded) 


Page  5 


Parameter 

TPR 

TT 

TTA-1 

TTB-2 

U 

U1 

VANEI 

VANE<|> 

WAE 

WAW 

WA/WT 

W8-1 

Wind-off 

W2 

W2R 

W2R-FS 


Tunnel  pressure  ratio 
Free-stream  total  temperature,  °R 
Primary  tunnel  total  temperature,  °R 
Back-up  tunnel  total  temperature,  °R 
Bypass  door  position 
Bypass  configuration 

Inboard  flow  throttling  vane  angle,  deg 
Outboard  flow  throttling  vane  angle,  deg 
East  wall  angle,  deg 
West  wall  angle,  deg 

Mass  flow  ratio  of  tunnel  make-up  air 
Wing  configuration 

Part  and  point  number  of  current  wind-off 

Engine-face  weight  flow,  lb/sec 

Engine- face  weight  flow  corrected  to  standard 
conditions,  lb/sec 

Full  scale  engine-face  weight  flow  corrected 
to  standard  conditions,  lb/sec 
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APPENDIX  C 

LIST  OF  COMMON  SYMBOLS  AND  CONVENTIONS 


d shock  tube  diameter 

IDA  average  distortion  index  . 

IDC  total  circumferential  distortion  index  I 

IDL  stall  margin  allocation  ratio  ' see  Tab^e  5,1 

IDR  total  radial  distortion  index  { 

IDT  total  distortion  index  ) 

M Mach  number 

PSi  overpressure  at  i'th  claw  probe  (i=l,2,3) 

Ap  shock  overpressure  at  blast  arrival  (scaled  to  1 atm.  ambient 
pressure  conditions) 
p pressure 

p pre-blast  tunnel  ambient  pressure 

p total  pressure 

Pt  '":e-blast  tunnel  total  pressure 

p 0 a’"-  age  engine  face  total  pressure  (psf) 

t2 

R2I  average  engine  face  total  pressure  for  inboard  inlet  (nondimensional) 
R20  average  engine  face  total  pressure  for  outboard  inlet  (nondimensional) 
r distance  from  end  of  shock  tube 

3M  fan  stall  margin  (Sec.  8.3) 

TT  pre-blast  tunnel  total  temperature  ( R) 

W2  engine-face  weight  flow 

W2R  engine-face  weight  flow  corrected  to  standard  conditions  (for 
full  scale  inlet  (lb/sec) 

W2R  = W2x(TT/519)  ' /(P  /2116) 

B side-slip  angle  2 
0 polar  angle  from  shock  tube  axis 

<*>  angle  between  shock  tube  axis  and  wind  tunnel  axis 
4>  blast  intercept  angle  (between  a normal  to  the  blast 
front  and  the  inlet  longitudinal  axis) 

Notes : 

In  most  piessure  plots  in  this  report  the  ordinate  label  designates 
the  variable  measured  and  the  vertical  scale  is  either  pressure/pt  or 

pressure/p^,  the  latter  being  designated  by  an  asterisk  after  the  ?abel. 

All  blast  pressure  values  in  this  report  designated  as  Ap  are  scaled 
to  a tunnel  ambient  pressure  of  one  atmosphere.  The  actual  test  values 
can  be  obtained  by  dividing  Ap  by  the  ambient  pressure  in  atmospheres 
(obtained  from  P in  Appendix  B) . 

Inlet  mass  flow  rates  in  this  report  are  always  expressed  in  terms 
of  full  scale  inlet  conditions.  To  obtain  model  values  divide  by  100. 

A separate  table  of  symbols  is  provided  for  Appendix  B in  that 
appendix. 
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